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Abstract 

High energy fixed angle scattering is studied in matrix string theory. The saddle point 
world sheet configurations, which give the dominant contributions to the string theory am- 
plitude, are taken as classical backgrounds in matrix string theory A one loop fiuctuation 
analysis about the classical background is performed. An exact treatment of the fermionic 
and bosonic zero modes is shown to lead to all of the expected structure of the scatter- 
ing amplitude. The ten-dimensional Lorentz invariant kinematical structure is obtained 
from the fermion zero modes, and the correct factor of the string coupling constant is 
obtained from the abelian gauge field zero modes. Up to a numerical factor we reproduce, 
from matrix string theory, the high energy limit of the tree level, four graviton scattering 
amplitude. 
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1. Introduction 



Weakly coupled perturbative string theory is conjectured to be described by the strong 
coupling limit of the matrix string super symmetric gauge theory [|l| [|| (see also Q for 
a review). In it is observed that the strong coupling limit corresponds to the infra red 
limit of the theory, and that symmetry more or less constrains the low energy effective 
action to be that of the light cone type IIA CFT. Matrix string theory, however, is a non- 
perturbative description of string theory and will contain phenomena, inaccessible to the 
infra red fixed point reasoning of , which can only be investigated by studying the gauge 
theory. A necessary first step in this direction is to derive, directly from the gauge theory, 
standard string perturbation theory. Although this program has not yet been realised 
several pieces of the puzzle have been found. 

Matrix string theory has classical solutions corresponding to interacting string world 
sheets in light-cone gauge |P(see also 0). The strings split and join via instanton like 
field configurations consisting of two regions; a core region, around the interaction point, 
where the fields do not commute, and an asymptotic region, away from the interaction 
point, where the fields do commute and where the eigenvalues glue together to form the 
Riemann surfaces of light-cone string theory 0(see also P]). Arguments have been given 
that in the large or the strong coupling limit [Q] the effective theory is indeed the 
type IIA CFT defined on the corresponding Riemann surface. However there are subtleties 
that mean that these arguments are not completely justified. An analysis of the effective 
expansion parameter for a perturbative gauge theory loop expansion [rU| shows that the 



loop expansion diverges for most physical scattering processes. There is thus, at present 
no rigorous approach to calculating the strongly coupled gauge theory in most situations 
of physical interest. 

A simple analysis does, however, lead to an explanation of the factors of the string 
coupling constant Qs associated with a string diagram 0. In the strong coupling limit 
there is, in addition to the Green-Schwarz light cone string action, a completely decoupled 
abelian gauge theory. The power of Qs is entirely due to the zero mode structure of this 
abelian gauge-field; it is just the difference between the number of closed, non-exact one- 
forms (the gauge field zero modes) and the number of ghost zero modes (exactly one, the 
constant mode). 

In this paper we focus on a scattering process where perturbative gauge theory cal- 
culations can be justified. The scattering process chosen is high energy, fixed angle, string 
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scattering, originally studied, from the string theory point of view, in [11||12] and re- 
analyzed from the perspective of matrix string theory in 0. Indeed it was originally 
pointed out in that this is a process where a perturbative gauge theory calculation 
might be justified. The scattering amplitude for high energy, fixed angle scattering is 
dominated by a classical world sheet configuration whose size, in the target space, is pro- 
portional to the incoming momenta. In other words high momenta correspond to large 
(in the target space sense) world sheets. By going to sufficiently high momenta, i.e. by 
making the world sheets sufficiently large, we can justify using a perturbative gauge theory 
calculation. This is similar to matrix theory [|13| graviton scattering calculations where 
one pulls the gravitons sufficiently far apart to be able to use a perturbative calculation. 

We will study explicitly the tree level contribution to the four ground state scatter- 
ing amplitude. Starting from the instanton-like classical matrix string configuration for 
such a world sheet we will estimate the fluctuation determinant and calculate exactly the 
contribution of the bosonic and fermionic zero modes. We will see that the zero mode 
calculation leads to all of the expected structure of the scattering amplitude. In addition 
to the correct power of Qs (observed in [0) we will reproduce the precise ten dimensional 
Lorentz invariant kinematic structure of the amplitude. The fluctuation determinant can 
not be calculated exactly, but we can nevertheless make some qualitative assessment of its 
overall form. The final result being that we will reproduce, up to an unknown numerical 
factor from the determinant, the four graviton scattering amplitude. 

Although we explicitly focus on the four graviton scattering amplitude the formalism 
is general and can be applied to the scattering of any ground state particles. The calcula- 
tion can also be generalized, in principle without difficulty, to loop scattering amplitudes, 
although we do not pursue this in this paper. 

We begin, in section 2, by covering some of the essential material necessary for the 
calculations of later sections. We start with a very brief review of matrix string theory. 
We then turn to saddle point world sheets in high energy string scattering and how they 
occur in light cone string theory. Next we review the instanton like field configurations 
corresponding to finite string interactions and how they embed into a solution correspond- 
ing to a string worldsheet splitting and joining at the interaction points. Finally we study 
the kinematics of four particle scattering in light cone gauge in the centre of mass frame 
and derive a few simple kinematic identities useful for section 8. 

In this paper we perform a one loop fiuctuation calculation around classical world 
sheet configurations of matrix string theory. Section 3 discusses the validity of such a 
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calculation and defines the Euclidean one loop action. In particular the Minkowski space 
Majorana-Weyl fermions are combined into complex fermionic coordinates and momenta 
which then allow a Wick rotation into Euclidean space. 

Section 4 analyses the fluctuation determinant. Sections 5 and 6 are devoted to the 
zero modes. In section 5 we focus on the behaviour of the zero modes in the neighbourhood 
of the instanton-like conflgurations, the objective being to understand how the core of the 
instanton, where the fields do not commute, effects the finiteness of the zero mode field 
configurations. Far from the core of the instanton these field configurations commute with 
each other and can be glued into global zero mode configurations for the global classical 
field configurations corresponding to light-cone string diagrams. These global zero modes 
we construct in section 6. 

In section 7 we discuss the incoming and outgoing wavefunctions and show how to con- 
struct ground state wavefunctions. In particular we explicitly construct from the fermionic 
coordinates of section 3 the graviton wavefunction. 

Section 8 is then devoted to the four graviton scattering amplitude, we explicitly inte- 
grate over all zero modes, the integration over the gauge and ghost zero modes leading to a 
factor of Qg (as argued by 0) and the integration over the fermion zero modes reproducing 
the ten dimensional Lorentz invariant kinematic factors for four graviton scattering. 

Having focused in previous sections on four graviton scattering we summarize in sec- 
tion 9 the general procedure for calculating high energy scattering amplitudes in matrix 
string theory, and point out the close connection with standard light cone gauge super- 
string calculations. We finish section 9 by discussing how the analysis might be performed 
away from the high energy limit, and speculate on how the calculations could be be given 
a rigorous basis. 

Technical details of the Euclideanization of the fermions and the fermion zero mode 
integrations are contained in the appendix. 

2. Matrix string theory, world sheets and high energy scattering 

2.1. Matrix String Theory 

Matrix string theory [000 is equivalent to ten dimensional SYM theory dimension- 
ally reduced to two dimensions : 

(2.1) 
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A careful derivation of this action (using the ideas of [n][15|[^) from the original matrix 
theory proposal |[T^] can be found in 0. All fields are N x N hermitean matrices. The 
index / for the bosonic fields runs from 1 to 8 and corresponds to the transverse directions 
of the ten dimensional target space. The 16 component fermion fields S split into S"" and 
S"", the 8s and 8c representations of SO{8). The matrices F'^, = 1, ■ ■ ■ , 9 are the spin{8) 
gamma matrices with in the Dirac operator p equal to the sixteen by sixteen unit 
matrix. Qs is the string coupling constant and the coordinate a runs from to 2n. The 



action (|2.1|) is conjectured to describe non-perturbative type IIA string theory compactified 



on a light-like circle with the number of quanta of momenta along the compactified 
direction. 

String world sheets are described in matrix string theory by commuting matrix con- 
figurations in which the A^ eigenvalues of the matrices form a branched covering of the 
cylinder, and hence form the Riemann surfaces of interacting light-cone string theory. 
These surfaces are characterized by strings of different length which split and join, the 
total length of the strings, which corresponds to the light-cone momentum, being pre- 
served. In the limit N oo the moduli space of the branched coverings is equivalent to 
the moduli space of all possible two-dimensional Riemann surfaces. A crucial ingredient 
of these configurations is that they are associated with a topologically non-trivial two di- 
mensional gauge field, which via Wilson lines generates the correct monodromy around 
all the branch points [H. Specifically a single valued matrix description of a multivalued 
branched covering of the cylinder is given by 

X = Udiag{xi,---,XN)U^, A^ = igMHdc.U), (2.2) 

where the eigenvalues Xi form the branched covering and the unitary matrix U (which is 
multivalued) generates the monodromies around the branch points. 

The gauge field configurations however are singular at the branch points and lead to a 
delta function singularity in the field strength. As was realised in this singularity can be 
resolved by instanton like field configurations. Far from the core of the instanton the fields 
commute with each other and their eigenvalues describe a Riemann surface with a branch 
point. Close to the core of the instanton, however, the fields no longer commute and the 
Riemann surface interpretation breaks down. It is to be expected that any classical string 
world sheet will have a corresponding classical solution in matrix string theory. 



Fig 1. Light cone string world sheets in matrix string theory 



In it was shown that there exists a physicaUy very interesting class of classical world 
sheets that preserve a certain amount of symmetry around the branch points and hence 
allow a simple construction of the instanton like field configurations. These world sheets 
are the saddle point classical world sheets dominating the amplitudes of high energy, fixed 
angle string scattering. Before reviewing the construction of the matrix string solutions 
we briefiy recall how these saddle point world sheets appear in light cone string theory. 



2.2. High energy scattering and light cone string theory 



High energy, fixed angle scattering processes were studied in string theory in [|TT| 
String theory simplifies enormously in this limit as the integral over world sheets localizes 
around a finite number of saddle-point configurations. It was conjectured in [|12| that 



the saddle point world sheets of different genus follow the same target space path, up to an 
overall scaling factor. This simplification permits one to study string perturbation theory 
out to arbitrarily high orders. One of the results of this analysis is that the perturbation 
series is extremely divergent, with the genus g world sheets contributing a factor of g^^ . A 
non-perturbative description such as matrix string theory should provide a natural cut-off 
to this divergence. First steps in an analysis along these lines were taken in . We will be 
less ambitious in this article and focus on retrieving from matrix string theory standard 
string perturbation theory. For the purposes of this article, we focus on high energy 
scattering as a means to justify using a perturbative gauge theory calculation. Below we 
describe how the saddle-point world sheets arise in light cone string theory. 



In light cone string theory the Virasoro constraints have been imposed at the classical 
level to eliminate the non-physical degrees of freedom. Specifically light cone string coor- 
dinates, = l/\/2{X^ ± X'^) have been eliminated from the dynamics by imposing the 
gauge 

X+ = r, (2.3) 

where r is the world sheet time direction. The coordinate X~ is then a function of the 
transverse coordinates X^ {I = I5 ■ ■ ■ 5 8) and the fermion fields. 

The light cone coordinates w = t + ia are defined using the Mandelstam mapping 
from the light cone diagram to the Riemann surface with uniformization z 



w 



Y^ptGiz.z.), (2.4) 



where G is an abelian differential with purely imaginary periods and real part equal to the 
Greens function. The existence and uniqueness of such a differential for arbitrary genus 



has been proved in |[T^] |jT^ . The branch points of the light cone string diagram are situated 
at the zeros of uj, i.e. the stationary points of w as a function of z. In the neighbourhood 
of a simple zero of u situated at zq we have w — wq = {z — zq)'^ . 
For tree level scattering we have 

w = T + ia = ^eiNilog{z - Zi), (2.5) 

i 

with z defined in the complex plane and e^A^i the p"*" momentum of the ith string. The 
are equal to +1 for incoming states and —1 for outgoing states. The lengths of the strings 
are proportional to their momenta. The light-cone interaction points are given by the 
zeros of dzW = 0, i.e. by the roots of the polynomial equation 

y = 0. (2.6) 

The classical solution for a string world sheet is given by 

X = ^eiPilog|z-Zi|, (2.7) 

where are the transverse momenta and z = z{w) is defined through ( p.5| ). 



For the classical field configurations (|2.7| ) the light cone action is given entirely by 
boundary terms 

S = -^j dMdcX'f = ^Y.', j daX'drX'l^^^, (2.8) 

i 

where the are the initial(final) times for the scattering process and the cr integral for the 
zth string runs over an interval of length A^^. For very big, i.e. far from the interaction 
region, equation ( p.5| ) can be inverted perturbatively and substituted into ( p.7| ) (|2.8| ) to 
obtain the results 

^ ^H^O [Pj - P^J^] log - 2, I + ■ ■ • 

' (2.9) 



2 ^ 



Zi Zj I ~\~ 



where the p momentum is given in terms of the transverse momenta pi and p"*" momentum 
iV^ by 

P7 = 1^ (2.10) 

The dots in ( p.9|) correspond to exponentially small corrections. The first term in the action 
( ^.9| ) is a phase factor corresponding to the time evolution of the incoming and outgoing 
states with respect to light cone Hamiltonian p~ , and cancels out in scattering amplitudes 
(see |]T^). The second term is the physically relevant Lorentz invariant contribution : 

^^ki.kjloglzi - Zj\, (2.11) 

where the ki are ten momenta defined to be incoming to the scattering process, i.e. ki = 
{Ni, \p(\'^ /{2Ni), pI) for incoming states and ki = —{Ni, \pl\'^ /{2Ni), p() for outgoing 
states. 

The saddle point contribution to the action is found by looking for the stationary 
point of the above term under variations with respect to the modular parameters Zi. 

2.3. Classical world sheets in matrix string theory 

The matrix string solutions for a general classical world sheet can be expected to 
be very complicated and to date no exact solution for a matrix string world sheet with 
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asymptotic states has been constructed ^. A successful strategy was however developed for 
the saddle point world sheet configurations which dominate high energy scattering ampli- 
tudes (see 0). Starting with the saddle point world sheet the authors of first identified 
the type of square root (holomorphic and/or anti-holomorphic) involved at the interac- 
tion points. They then constructed instanton-like field configurations whose behaviour far 
from the core matched the square root behaviour of the world sheet. The approach is an 
approximation, in that the instanton-like field configurations have a finite size and should 
strictly be matched onto more than just the dominant square root dependence on the world 
sheet coordinates. The true instanton-like field configuration would thus be expected to 
be slightly modified away the core. 

In this program was implemented explicitly for the construction of the four string 
scattering process. Below we briefly review the main points. In the centre of mass frame 
the scattering process describes a two dimensional plane and the corresponding classical 
world sheet can be described by just two of the eight transverse coordinates^, X'^ and X^. 
For convenience these are combined into a pair of complex fields defined by 

X = ^{X^ + iX^) and X* = ^{X^ -iX^). (2.12) 
V 2 V 2 

The light cone world sheet has two branch points Zq and Zq given by the two roots of 
the quadratic equation for four particle scattering of (|2.6|). The dominant classical world 
sheet has the property that the complex scalar field X is an anti-holomorphic, respectively, 
holomorphic function of z in the neighbourhood of the two branch points. Specifically : 

i 

and 

9^X\.-.-=y2^^L-.~=0- (2-14) 

i * 

This symmetry allows a simple construction for the instanton-like field configurations 
around the branch points. The starting point is the four dimensional self dual YM equation 

^ Exact matrix string world sheets have been constructed and studied in depth in ||6| [^] , however 
these are for classical world sheets which are entirely holomorphic or entirely anti-holomorphic 
functions of the world sheet coordinate and thus cannot represent the matrix string versions of 
the classical world sheets ( |2.7| ). We thus do not consider them further in this article. 

^ Scattering processes involving more than four particles would require three or more transverse 
coordinates to describe their classical world sheets. 
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dimensionally reduced to two dimensions. The two possible signs for the self dual equations 
correspond to the two different holomorphic behaviours of equation (|2.13|) . The instanton 



solution around the branch point Zq is obtained from F^i, = t^^paFpu with eoi29 = +1- 
In terms of the two dimensional fields the equation reads 

Fww = — 

D^X=0 ' (2-15) 
Di,X* =0 

where, in addition to the complex scalar fields defined in ( |2.12| ) we have introduced complex 
two dimensional gauge fields and coordinates 

A^^(A'^-iA^) and w ^ ^(a^^ + ia^). (2.16) 
\/2 V2 

A simple solution to these equations of motion corresponding to a simple branch point 
involving two eigenvalues is 

X = UXU^ and A = ig,U^dU with X = BViLt^, and t/ = e^^"^^^ (2.17) 

As already described this leads to a delta function singularity in the field strength at the 
interaction points. This singularity can be removed once we are working with complex 
coordinates X, by using a complexified "gauge" transformation G which also has a singu- 
larity at the origin tuned in such a way as to leave a singularity free field strength. With 
this insight the solution can be written in the form 

X =UGXG~^U^ 

(2.18) 

A = -ig,[G-\d^G) + U\d^U)], 



where the diagonal matrix X and unitary matrix U are as in (|2.17|) and the matrix G is 
given by 

Q^^a{ww)T,_ ^2.19) 

This ansatz automatically satisfies the last two equations of ( p.l5| ) with the first equation 
leading to a differential equation for a 

(a2 + -a,)a = -^rsmh2a with (2.20) 



where r = \/ww is the radial distance from the branch point. The boundary conditions 
are necessary for a finite solution. In particular the second boundary condition ensures 
that there are no ^ pole terms in the gauge field A and hence no delta function singularity 
in the field strength Fy^^j. 

There is no explicit solution to this equation but it can easily be solved numerically. 
In terms of the function a the expressions for the scalar fields X and the field strength 
Fuiiju are given by the simple expressions 

X = i?v^(coshQ! Ts + zsinhar2) and F^^, = rsinh2Q!. (2-21) 

9s 

For simplicity we do not include the final gauge transformation U 

The instanton like solutions ( |2.21| ) are embedded into a global solution, defined on 
the cylinder, corresponding to the classical world sheet. The parameter B of equations 
( p.l7 ),( p.20|) and ( 2.21 ) is determined in terms the four external momenta by 



1^,2 _ \PlPl-PlP3\{Ni+N2) 

' ' ^yNlN2N3N^ ' ^ ■ ^ 

where the pi and p* are the complex transverse momenta associated with the field X and 
its complex conjugate. The A^^ are the p"*" momenta for the ith string. 

Finally we note that there are two physical scales associated with the instanton. Firstly 
the differential equation ( |2.2C1| ) can be given a dimensionless form by absorbing the coupling 
constants into a rescaling of the radial coordinate. In other words the instanton has a 
natural world sheet scale : 

knst ~ ( ^ ) ■ (2.23) 

Secondly, a simple analysis of the limiting behaviour of the field X at the origin shows 
that there is a minimal target space "distance" between the two strings : 



d,nst = VT^r{X{0)X*{0))r^gi\B\-^. (2.24) 

2.4- Kinematic relations for four string scattering 

We will explicitly study the four string ground state scattering process. The calcula- 
tion is performed in the centre of mass frame. In this subsection we give various kinematic 
identities useful for the calculations of later sections. 

All particles are massless. They have transverse momentum lying in the plane 
specified by the complex numbers pi = 1/V2{p} + ipf). The p'^ momentum is given by 
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the length A^^ of the string. To make contact with standard conventions for scattering 
processes we define aU ten-vectors of momenta ki to be incoming momenta : 

, (2.25, 



where 



(-iV„-^,-p„-p*) for z = 3,4 



\P^\'=P^P: = liip}? + {p^r)■ (2.26) 



The dot product of two momenta then reads 

ki.kj = -pfpj - p'j'p~ + {piP* + p*Pj). (2.27) 
In the centre of mass frame conservation of momentum reads 

A^i + A^s = A^3 + = (2.28) 

where in the final line (conservation of p~ momentum) we have used the first and second 
lines (conservation of transverse and p"*" momenta) to simplify the result and to define the 
quantity p^. Inverting the third line we can write and |p3p in terms of p^ 

\Pi\ = -J^P \p''\ = -j^P ■ (2.29) 

Finally it is useful to define the ten dimensional Lorentz invariant quantities for the scat- 
tering process 

s = -{ki + k2f = -2/ci./c2 

t = -{k2 + ksf = ~2k2.k3 (2.30) 

u = -{ki + kaf = -2/ci./c3, 

where the sign conventions are, as before, for /ci,/c2,^3 and A;4 incoming momenta. Using 
the definitions for the ki ( p.25| ) along with the identities ( |2.28| ) and (p.29| ) the Lorentz 
invariants s,t and u can be written as 

s = 2p^ 

(N1N3 + N2N4) 2 2 
^ = -2 P -2q ^ith =pIp3+P3PI- (2.31) 

u = -2 — p +2q 
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in the above expressions is given in (|2.28| ) ( pT29|) . It is easy to check that the above 
expressions satisfy the identity relating s, t and u for massless particles 

s + t + ti = 0. (2.32) 



3. Perturbative calculations and the Euclidean action 

There are two standard limits one can take in matrix string theory to recover pertur- 
bative string theory. Firstly one can take the limit N ^ oo followed by Qs —>■ 0. By the 
matrix string theory conjecture this is the limit necessary to recover perturbative type IIA 
string theory in uncompactified ten dimensional space. Secondly one can hold finite 
and send Qs 0. By Susskind's conjecture this limit corresponds to string theory 



compactified on a light-like circle with the number of quanta around the compact 
direction. 

Since both are strong coupling limits one would not expect perturbative gauge theory 
calculations to be justified. Indeed, a perturbative gauge theory expansion would involve, 
through the three and four point vertices, arbitrary powers ofl/gg which diverge as Qs — > 0. 
However, when deciding whether or not a perturbative calculation is justified it is the 
effective parameter weighting the loop expansion that is important not the size of the 
coupling constants. This is determined by both the coupling constants and the masses 
of the particles exchanged in the propagators. The masses of the quantum fiuctuations 
are determined by the commutator term in the action ( |2.1[ ) and are hence proportional to 
I/qs- The mass of the propagators can thus compensate for the strength of the coupling 
constants. A good analogy to bare in mind is the top quark and Higgs in the standard 
model. The coupling of the top quark to the Higgs field is enormous (it is this that gives 
it its enormous mass). This does not mean that perturbative standard model calculations 
involving a virtual top quark and Higgs are not justified (ask any phenomenologist). There 
is precise cancellation between the mass in the propagator of the virtual top quark and 
the coupling of that virtual top quark to the Higgs. An identical cancellation happens in 
matrix (string) theory. 

In the context of matrix string theory an analysis of the balance between these two 
effects, in the large A^ limit, was carried out in |jT^. Starting from a bosonic background 
field there is a systematic expansion in powers of derivatives of the background fields for 
the effective action that one can calculate. It's overall form is entirely determined by 
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dimensional reasoning and the identification of tlie loop counting parameter [ 21 1 . The 
result is that the effective action can be written in the form 



5'eff = Ida 



F' + J2^L with /:L = E^^"' x4n+2L-4 - (3-1) 

n=2 



L=l 



where F^'^/X^"^ means bosonic terms with 2n derivatives in the numerator and 2m powers 
of the scalar fields in the denominator. 

The analysis of [TO] showed that the two standard limits mentioned above lead to 
divergent loop expansions. For the case N ^ oo before Qs ^ the tensor structures of 
the terms F2n^j^4n+2L-4 g^j^g such that the loop expansion diverges with positive powers 
of A^. The physical cause of the divergence is that neighbouring strips (eigenvalues) of the 
background long string configuration come arbitrarily close together in the limit N oo. 
The off diagonal field variables connecting the neighbouring strips thus become massless 
in this limit. Holding N finite and sending Qs leads to a different kind of problem. 
Naively from (^^) the effective action would be well behaved in this limit. However close 
to the interaction points the instanton-like field configurations depend upon Qg. Taking 
this into account leads to a loop expansion weighted by inverse powers of Qs, (see [|T0| ) 
which diverges in the limit Qs ^ 0. Physically the reason for this divergence is that, as 
can be seen from ( pj.'24D , the "minimal distance" between the two strings tends to zero in 
the limit Qs ^ 0. This again leads to massless fluctuations and a divergence. To be able 
to calculate in these two limits would thus first require the development of direct strong 
coupling techniques. 

In this article we will avoid these problems by studying a third limit which will allow us 
to calculate a perturbative string theory scattering amplitude by performing a perturbative 
gauge theory calculation. The limit corresponds to high energy string scattering in the 
finite version of matrix string theory. Specifically we will hold fixed, hold Qs fixed 
and study the dominant contribution in the limit in which the external momenta, pi, of 
the scattering process tend to infinity. Note that can be large and Qs can be small. 
What is important is that we send pi oo before taking any other limit. 

The strategy of this paper is to perform a fiuctuation calculation about a classical 
matrix string theory background corresponding to one of the dominant saddle point world 
sheets found by Gross and Mende |]11| [0 in the pi oo limit. As can be seen directly 
from ( p77| ) the target space size of the saddle point backgrounds are proportional to the 
external momenta, pi ^ oo is thus a limit in which the target space size of the background 
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becomes infinite. Looking at loop expansion (|3.1| ) we see that there are always more 
powers of the background field in the numerator than in the denominator. In the limit 
in which the background becomes infinitely large these terms will thus be scaled away 
altogether. Physically the limit has the effect of separating, in target space, the individual 
strips (eigenvalues) of the matrix background, with the result that the off-diagonal elements 
connecting together the strips become infinitely massive. Again this simple argument 
could be invalidated by the existence of the interaction points where two eigenvalues are 
connected by a branch point and come close together. As already stated the instanton- 
like field configurations lead to their being a "minimal distance" ( |2.24| ) between the two 
eigenvalues. Reading off directly from (|2.24|) ( |2.22| ) we see that this minimal distance tends 



to infinity if we take the limit pi ^ oo while holding and Qs fixed. There is thus no 
problem of massless fiuctuations as there was for the (7^ — > limit. 

To summarize, taking the limit — 00, before taking any other limit will lead to a 
loop expansion in which the contributions from higher order loops will be scaled away. It 
is thus justified in this limit to use a one loop calculation. As will be discussed below the 
calculation reduces to the evaluation of the determinant for the quadratic fluctuations (by 
the above arguments this will be effectively equal to one) and to an integration over the 
zero modes. 

A final comment is in order. Asymptotically far from the interaction points the clas- 
sical background splits into separate blocks for the different strings. Each block is propor- 
tional to a unit matrix. The off diagonal elements within a block are thus massless. This 
is a result of the fact that we are not using true matrix string wavefunctions. Presumably 
there is some LSZ type reasoning that could be developed to justify replacing the true 
matrix string wavefunctions by simple diagonal blocks. 

3.1. Euclidean fermions 

The fermions of the action ( |2.1[ ) are Majorana Weyl fermions in ten dimensional 
Minkowski space and consist of sixteen real fermionic components. We will need to work 
in ten dimensional Euclidean space (dimensionally reduced to two dimensions) where the 
instanton configuration and functional integral are defined, and where it is not possible 
to define Majorana Weyl spinors. In addition the Majorana Weyl fermions satisfy the 
commutation relations 

{S''{a),S\a')] = Tx5''^5{a-a') and {S''{a),S\a')] = T:5^h{a-a'), (3.2) 
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where a,b and d, 6 denote respectively the Sg and 8c representations of SO{8). In other 
words the spinor S"- is simultaneously a coordinate and its conjugate momenta, and sim- 
ilarly for S°'. The spinor variables can however be combined into complex spinor coordi- 
nates, which are distinct from their conjugate momenta, and which can be used to define a 
Euclidean action. The complexification procedure means that the SO{8) symmetry of the 
theory will no longer be manifest, although, as we will see, it is restored in the calculation 
of physical amplitudes. A similar complexification procedure has to be carried out in light 
cone superstring theory |T^ . 

The calculation of this paper is up to quadratic order only in the background fluctu- 
ations. To quadratic order the Dirac operator involves only four gamma matrices, and 
r^, which couple to the background gauge field and and coupling to the background 
fields for and X^. This fact permits a simple definition for the Euclidean action. We 
start by choosing a basis for the gamma matrices such that the four gamma matrices, 
F'^,F^,F^ and F'^, can be written as four by four blocks (see appendix). We then define 

four spinor coordinates 9^,6^,0^, 0^ and four corresponding momenta Xa,X^ 
(see appendix) 
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(3.3) 

The indices A and A take the values A, A = 1,2. Note that 9 is not the complex conjugate 
of 9, they are distinct variables. The notation is chosen to match standard conventions 
for spinors in four dimensional space (see [|2^). This choice of coordinates breaks the 
spin{8) symmetry of the spinors down to U{1) ® U{1) ® SU{2) ® SU{2). The two t/(l)s 
correspond to SO{2) rotations in the and (with the basis of gamma matrices chosen 

in the appendix) X^^X^ planes. The SU{2) ® SU{2) corresponds to an 5'0(4) for the 
remaining four dimensional space X^ ,X^,X'^ and X^ . The charges of the fermions under 
the rotations in the X^^X"^ and X^,X^ planes as well as the transformations under rotations 
in the X^, X^, X^, X^ can be read off directly from the form of the eight dimensional F-'^'^ 
(see appendix). The final result of this decomposition is that we have the following bosonic 
and fermionic coordinates. 



X,X' 



X,X' 



1 



-^(X^±iX^] 
V2 



^2' 2' 



V2 



{X^ ± iX' 
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(3.4) 
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The first and second arguments of tfie fermions are tlieir ^7(1) cfiarges under rotations in 
respectively tfie X'^,X'^ and planes. The absence or presence of the bar above 

the fermions indicates how they transform under rotations in the , , X''' , space. 
Specifically under such an infinitesimal rotation we have 

66^ =\w^-{a^-)^Be^ with (a-)^ = ^(a-a"-a-0 (3.5) 

The indices mnpq run from 5 to 8. The a and a matrices are given in terms of the pauli 
matrices by 

(a-)^^ = (z,ri,r2,r3) and (a-)^^ = (-z, n, r2, rg). (3.6) 

Note that this decomposition of the fermions into coordinates and momenta differs from 
the standard decomposition used in light cone superstring theory [|19[ where one breaks 



the spin{8) symmetry of the spinors into U{1) ® SU{4). The classical background further 
breaks the SU{4) down to U{1) x SU{2) x SU{2). 

In section 8 we will use the transformation properties of the spinors 9,9,9 and 9 
to construct incoming and outgoing wavefunctions transforming under SO [8) rotations. 
Specifically we will construct the combinations of 9 and 9 corresponding to the X ,X* ,X ,X* 
and X'^ components of an SO (8) vector and similarly for 9 and 9. These "left" and "right" 
SO [8) vectors will form incoming and outgoing graviton states. We will then explicitly 
calculate the four graviton scattering amplitude. 

3.2. The Euclidean action 

It is convenient in this section to use ten dimensional notation for the bosonic fields. 
Indices fi run from |U = 0, ■ ■ ■ , 9 and split into indices 0, 9 for the two dimensional cylindrical 
coordinates on which the fields depend and / = 1, ■ ■ ■ , 8. All bosonic fields are denoted by 
with A^ = X^ . We then split the fields into a background part, A, corresponding to 
the classical matrix string solution and a fluctuation part V : 

Atotai = A + V, (3.7) 

with = . To calculate the effect of quantum fluctuations about the classical con- 
figuration the action needs to be gauge fixed. We use the standard background covariant 
gauge fixing term 

Cgf = iD^V^f with D^ = d^ + -[A^,]. (3.8) 



9. 



s 
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The Euclidean Lagrangian for tlie quadratic fluctuations reads 



Ceu = Tr [-(L>^y")2 - 2 — V^[F^\ F"] + X^pO + c'D'^c] . 

9s 



(3.9) 



where V are the bosonic fluctuations and 9 and A the fermionic coordinates and momenta 
defined in ( |3.3| ). The fields c are the ghosts. D is the background covariant derivative 
given in ( |3.8| ) . Using the definition of gamma matrices given in the appendix the fermionic 
part of the Lagrangian reads 



( 

-MX, 



[X' 



-Mx, 

9s I- ' 



\ 



(3.10) 



4. The fluctuation determinant 

The fermionic coordinates A act as Lagrange multipliers and integrating over them 
leads to 5{p9). In other words the functional integral for the fermionic coordinates 6 is 
projected down onto their zero modes. This is discussed in more detail in the context of 
light cone string theory in [25|. The Jacobian factor from this projection is (iet(lf)^If)^ . 

In this section we focus on the contribution from the non zero modes i.e. on the 
determinants. For the four string scattering background considered in this paper both the 
fermion term and the boson term consist of four by four non-trivial blocks. Canceling 
part of the trivial block of the bosonic determinant with the ghost determinant the total 
fluctuation determinant J reads 



J 



(4.1) 



I- ' J/4x4 

where the indices p, a take the values 0, 9, 1, 2, and zero modes have been excluded from 
the determinants. The fermion determinant is in four by four spinor space and the boson 
determinant in four by four p, a space. 

It is well known that the fermionic and bosonic fluctuation determinants precisely 
cancel in a self dual background. Using the definition of gamma matrices given in the 
appendix it is only a short calculation to verify that this is indeed the case. Specifically 
for the instanton background of section 2.2 we have 
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(4.2) 
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where we have we have used complex indices, i.e. rows labeled from top to bottom by 
w,x,w,x and columns labeled from left to right by w,x,w,x. Using ( p.lO|) we have 



/ ^ 

p^p = hD^ + -[F, ]) with F 
^ 9s 



^ WW J-^tn-^ 



(4.3) 



V -DuqX F^ww / 



leading to J = 1. 

In other words if the background was self dual everywhere the fluctuation determi- 
nant would be equal to one. The matrix string background corresponding to four string 
scattering is however only locally self dual around the branch points and its fluctuation 
determinant would be corrected by perturbations from the non self dual part away from 
the interaction points. The important point to note is that, in the regime where it is 
justified to use a one loop calculation, these are just small corrections and to lowest order 
can be dropped. The conclusion from this analysis is that the interaction points do not 
lead to a singular contribution to the determinant. This is true even in the limit where Qs 
and hence (via (|2.23| )) the size of the instanton, tend to zero. 

It is interesting to compare this result with what one would expect from a CFT 
fiuctuation determinant around a branch point. In light cone super string theory 6 
and A have the conformal weights 1 and 0, respectively. Potentially there is a singular 
contribution coming from the square root cut point which hides a world sheet curvature 
singularity. Smoothing this out over some cut off distance e leads to a singular contribution 
l/eT^ coming from the induced Liouville action for the bosons and fermions. However for 
the choice of weights given above the bosonic and fermionic central charges cancel and 
there is no singular contribution. This can also be seen directly from a comparism of the 
determinants (see P^). 

The conclusion of this section is that there are no singular contributions to the deter- 
minant coming from the interaction points and that this is consistent with the light cone 
superstring calculation. This does not mean, however, that the fiuctuation determinant 
will be equal to precisely one. We will return to this point in section 8. The remaining part 
of the functional integral is an integration over a finite number of zero modes and collective 
coordinates of the background configuration. In the next two sections we construct these 
modes. 
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5. Instanton zero modes, local considerations 

In this section we focus on the field configurations in the neighbourhood of the in- 
stanton. The objective is to understand how the core of the instanton, where the fields 
do not commute, effects the finiteness of the zero mode field configurations. Far from the 
core of the instanton the fields commute and the zero modes correspond to those defined 
on a Riemann surface in the neighbourhood of a branch point. In particular we will see 
that there are zero modes which, from the Riemann surface point of view, would appear 
to diverge at the branch point, but which are rendered finite by the non-commuting core 
of the instanton. 



These "singular" modes will play a crucial role in the next section where we construct 
global zero modes. The instanton is embedded into a global classical solution corresponding 
to a classical string world sheet. The instanton zero modes will thus be embedded into 
global zero modes which tend, asymptotically far down the strings, to constant values. 
From the point of view of the world sheet these modes can be written as holomorphic 
and/or anti-holomorphic functions of the cylindrical coordinates. The singularities at the 
branch points are important since it is only with singularities that one can construct non- 
trivial global modes. 

The global zero modes will include the moduli for the classical string world sheet along 
with their superpartners and zero modes for the abelian gauge theory defined on the string 
world sheet. 

For the purposes of this section we will treat the instanton as being embedded into the 
complex plane. Below we will explicitly construct the bosonic and fermionic zero modes, 
imposing the condition that they must be finite everywhere (except at infinity). 



1 



1 




Fig 2. Finite field configurations for "singular" zero modes 
shown in the w plane and the z = y/w plane. 
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5.1. Bosonic zero modes and collective coordinates 

The classical solution for the instanton (|2.18| ) ( |2.19| ) ( p.20|) involves only the two dimen- 
sional gauge field A and two of the 8 possible scalar fields, Xi and X2. There are three 
types of bosonic zero modes for the fiuctuations about this configuration characterized by 
which fields they tend to asymptotically. Firstly there are the two-dimensional gauge field 
zero modes. These come from the fact that the gauge fixing term doesn't completely fix 
the gauge. In other words there are local gauge transformations that are zero modes of 
the gauge fixing term. The second correspond to translations and deformations of the 
instanton solution, and asymptotically involve only the scalar fields Xi and X2. Finally 
there are modes involving X3, ■ ■ ■ , Xg. All three types are zero modes of the eigenvalue 
equation for quadratic fluctuations about the classical instanton background : 

[D,D,rj^, + ^[F^., = A^-VJ'^). (5.1) 

In this equation the vj^^ are the fluctuation modes and A*^"^^ their eigenvalues. For com- 
pactness the equation has been written in its ten dimensional form. The ten dimensional 
indices fx.jy = 0, ■ ■ ■ , 9 split into the two dimensional indices a, (3 = 0, 9, corresponding 
to the cylindrical coordinates, and the indices / = l,---,8 corresponding to the scalar 
matrix flelds. The covariant derivatives Df^ are covariant with respect to the background 
flelds, i.e. -D^ = + i/g[A^, ], and F^, is the background fleld strength. The standard 
Feynman background gauge fixing term {Dfj,A^)^ has been used. 

5.2. Gauge field zero modes 

For the gauge field zero modes we search for gauge-fields, V^, that are locally pure 
gauge, 

= D^A, (5.2) 

and which satisfy equation (|0|). Substituting ( p. 21 ) into ( p.l| ) and commuting the covariant 
derivative to the right through the term cancels the F^i, term and leads us to look for 
modes A that satisfy 

DpDpA = 0, (5.3) 

i.e. to look for zero modes of the gauge fixing term. The resulting gauge fields must be 
finite everywhere (except possibly at infinity) and must tend to some function times ts at 
infinity so that they commute with the bosonic fields. Note that equation(|5.3|) is just the 
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equation of motion for the complex bosonic field X. So we can immediately write down a 
candidate zero mode : 

A(^)=X*. (5.4) 
with X given by ( p.21| ). The gauge field zero mode l/^"^) = D^A*-^-* is given by 

— — _ |^(-cosha + rc^rCK sinha)r3 + z(-sinha + rdrOi cosh.a)T2\ ■ (5-5) 

One still has to check that l/^"^) is finite at the origin, since it potentially diverges as 1/r. 
Using the fact that for small r a{r) = — l/21ogr + a + (9(r^), we see that it is indeed finite 
at the origin : 

ri 1 

y^"^) ^ ^-e""[r3 + zr2] as r ^ oo, (5.6) 
u)2 2 

and is thus a valid zero mode. 

The most general solution for the gauge zero modes can be found by writing equa- 
tion(5.3) in terms of the two-dimensional fields and covariant derivatives : 



DpDpA = D^D^A + D^D^A - ^{[X, [X*,A]] + [X*, [X,A]]) = 0. (5.7) 
It is straightforward using ( |2.15| ) to see that it is satisfied by 

A(-) = a(") + (A("))*, (5.8) 

with 

^{n) ^ i w^'-^X* ^ w'^Ts for n = i,f,f,--- 

\«;"l2x2 for n = 1,2,3, ■■■ ' ^ ' ^ 

and 

vW=L»^A("-^). (5.10) 

All these solutions are of course written in the singular gauge in which asymptotically far 
from the instanton all fields are diagonal. 

The physical meaning of these zero modes is given by their behaviour far from the 
non-commuting core of the instanton where they can be interpreted as zero modes of an 
abelian gauge field defined on a Riemann surface in the neighbourhood of a branch point. 
They are zero modes of the abelian gauge-fixing term (9q,Vq,)^. We see that we have the 
set of abelian pure gauge zero modes. 
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Vw=w ,l,w^ ,w,w^ jw"^ , ■ ■ ■ , (5.11) 
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in the complex w plane. In particular there is a mode which is finite at the origin 

of the instanton but asymptotically behaves as w~^. Unwinding the branch cut w = z^, 
Vz = 2y/wVw we see that we have the full set of abelian gauge field zero modes, 

Vz = l,z,z^z\z\---, (5.12) 

defined in the complex z plane. In section 4 we study the zero modes about the classical 
solutions corresponding to light-cone Riemann surfaces. In other words we embed the zero 
modes constructed above into global abelian zero modes defined on the light cone Riemann 
surfaces. The fact that we have the full set of zero modes ( p. 121 ) means that we will be 



able to construct arbitrary abelian zero modes on the Riemann surface which are not in 
any way pinned to, or constrained by, the interaction points. 

The conclusion from this section is that for the construction of the global zero modes 
for the complex gauge field V we are allowed a, {w — Wq)~^ ~ l/{z — Zq) singularity at 
z = Zq and a {w — Wq)~^ 1/{z — Zq) singularity at z = Zq . 

5.3. Zero modes for the ghosts and for X^, ■ ■ ■ , 

The zero mode equation for the ghosts is identical to equations (p.3|)( |5.7D . This allows 
to automatically write down all their zero modes. They are given by the solutions for 
A (equations (|5.8|) (|5T9|)) of the previous section. The important point is that there is no 
mode that asymptotically behaves as and is finite at the origin, so there will be no 
non-trivial global zero mode we can construct. 

Since the background only depends on the scalar fields X^ and X"^ the zero mode 
equation for the scalar fields X^, ■ ■ ■ ,X^, contains no F^i, term and is thus also identical 
to equations (|5 . tj|) (p77|) . So again there will be no non-trivial zero modes. 

5.4. Zero modes for X^ and X"^ 

The most obvious zero mode to write down for the fields X = l/^/2{X^ + iX"^) and 
X* is that corresponding to the translation of the instanton configuration. To construct 
this mode one proceeds as for the four dimensional instanton by taking the derivative of 
all the fields and at the same time shifting them by a gauge transformation so that the 
background gauge fixing condition, D^V^, is preserved. The result being that the zero 
modes of ( |5.1D that correspond to translations are given by 

= d^A^ - D^A^ = F^^. (5.13) 
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The first term corresponds to translation in the a direction, the second term is the gauge 
transformation with gauge parameter A^. Substituting ( |5.13|) into (|5.1J ) and using the 
equations of motion and the Bianchi identity this is indeed seen to be a zero mode of ( |5.1| ). 
Written in terms of the two dimensional complex fields of section 2 the solution reads 

Its asymptotic behaviour, indicated on the r.h.s, shows that from the world sheet point 
of view it is singular at the origin. The non-commutative core of the instanton, however, 
renders the field configuration finite at the origin (see equation (|5.6|) ). As already stressed 
the fact that it is a singular zero mode from the world sheet point of view will allow us 
to construct a non-trivial global zero mode. It will correspond to deforming the classical 
world sheet by moving a branch point. 

It is important to check to see if their are any other singular zero modes. The only 
other possibility would be a mode asymptotically behaving as w~^. We write equation 
( ^.1|) in terms of the two dimensional covariant derivatives, field strengths and scalar fields : 




+ 2- ( \ ^^,1 ' J)^ (;)=0, (5.15) 



where DpDp, in terms of two dimensional covariant derivatives and fields, has already been 
given in ( p.7| ). It is possible to find simple generalizations of ( |5.14| ), satisfying ( |5.15| ), for 
arbitrary half-integer powers of lu : 

along with the trivial solution 

Y = iD"l2x2 n>0. (5.17) 

However there is no obvious, simple generalization of ( |5.16| ) for powers o( w. It would 
seem natural, however, that there would NOT be a mode behaving as which is also 
finite at the origin. 

We will thus assume that this is the case (to prove it would mean analyzing carefully 
the coupled differential equations of ( |5.15| )). With this assumption the conclusion from 
this section then is that for the construction of the global zero modes we are allowed a 
{w — Wq)~^ ~ 1/{z — Zq) singularity aX z = Zq and a {w — Wq)~^ ~ 1/{z — Zq) singularity 
at z ^ Zq . 
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5.5. Fermionic zero modes 

The fermion zero mode equation reads 

lpe = 0, (5.18) 

with^^ given explicitly in ( |3.1C1| ). The most important zero mode is that coming from the 
Euclidean version of the infinitesimal Minkowski space supersymmetry transformation, 

= F^.r^.e, (5.19) 

where e is a constant ten-dimensional Majorana-Weyl spinor. This can be show to be a 
solution of (|5.18| ) by using the gamma matrix identity TpT^,^ = Fp^,^ + l/2(?7p^r,^ —r]pt,r^), 
followed by the Bianchi identity and the equation of motion for F^^. 

It is easy to see from ( p.lO|) that the Euclidean fermion zero mode equivalent to ( |5.19| ) 
has the form 9 = 9 = and 

where 9^ and 9 are two component constant, complex spinor variables. The asymptotic 



behaviour on the r.h.s. of ( |5.20| ) shows that from the world sheet point of view the mode 



is singular at the origin. As for the bosonic zero modes the field configuration is in fact 
finite at the origin. 

Again it is important to check that there are no other "singular" zero modes. Using 
the explicit expression for p of ( |3.10| ) it is relatively simple to see that 9 and 9 have the 
same solutions, respectively, as A*^") and A*^") (equations (|5.7|) (|5.8|)). In other words there 
is no zero mode for 9 or 9 that from the world sheet point of view appears to have a 
singularity at the branch point but is nevertheless finite at the origin. 

For 9 and 9 we have already found the "singular" zero mode ( |5.19| ) and there are no 
other possibilities (we cannot, for example, have 9 behaving as since asymptotically 
9 must be an anti-holomorphic function of w). 

The fermion zero mode ( p.20|) is for the self-dual configuration around the branch 



(5.21) 



point w = Wq . The zero mode around the anti-self-dual configuration'^ of the Wq branch 
point can also be read off from ( |3.1C1| ). It is given by ^ = 6* = and 

9\^fF^^ D^X*\f9P\ f9Pw-^ 
9 J \D^X F^^ J {d")^ \9Pw-^ 

The conclusion is that for the construction of global fermion zero modes in the next 
section 9,9,9 and 9 can respectively have the singularities 1/{z — Zq), 1/{z — Zq), 1/{z — Zq), 
and l/{z — Zq). 



For the anti-self-dual instanton we have D^X* = 0, D^X = and Fww + -^[X* ,X] 
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6. Zero modes, global considerations 

In the previous section we analyzed in detail the bosonic and fermionic zero modes in 
the neighbourhood of the instanton configuration for the branch points of the matrix string 
world sheet. In particular we identified zero mode field configurations which are finite at 
the branch points but from the world sheet point of view appear to have a singularity at 
the branch point. In this section we take these zero modes and embed them into non- 
trivial zero modes for the gauge, boson and fermion fields for global classical solutions 
corresponding to light cone Riemann surfaces. 

We will construct global zero modes which tend to constant values asymptotically far 
down the strings. Specifically we will construct modes satisfying the boundary conditions 

Vr = drV^ = 0, drY^ = 0, = for r = Ti or Tf, (6.1) 

where ri and Tf are the initial and final times. 

From the world sheet point of view all zero modes can be decomposed as holomorphic 
and/or anti-holomorphic functions of z, with possible simple pole singularities at the branch 
points. At the end of sections 4.2,4.3,4.4 and 4.5 we enumerated the possible pole structures 
for the bosonic and fermionic modes. In the following sections we use these singularities to 
explicitly construct non-trivial global zero modes satisfying the boundary conditions ( |6.1| ). 

6.1. Abelian gauge and ghost field zero modes 

The two dimensional gauge field V = 1/V2{Vt- — iVa) must be a holomorphic function 
of w or equivalent ly of z. We have already seen in the previous section that it can have 
the poles l/{z — Zq) at the branch point Zq and l/{z — Zq) at the branch point Zq . This 
allows us to construct the following zero mode 

E^e^I^^o^ for z^z,, (6.2) 
with the Gi real numbers satisfying 

$^Qa, = 0. (6.3) 

i 

Note from equation ( p.6|) that the numerator of 
hence gives the correct pole structure 

Asymptotically far down a string {z Zi) the boundary conditions ( |0|) are satisfied 
with Vcr tending to the constant value ai/Ni as indicated in ( |6.2| ). A Wilson loop encircling 
the string thus generates the U{1) element es^"'. 

The mode ( |6.2| ) is the most general abelian differential defined on the string diagram 
Riemann surface. 

For the ghost field zero modes no singularities are allowed at the branch points and 
there is thus just a single, globally constant zero mode. 



6.2) has zeros at the branch points and 
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6.2. Bosonic field zero modes 

The non-trivial bosonic zero modes satisfying the boundary conditions (|6.1| ) involve 
only the bosonic field X = l/\^{Xi +1X2) and it's complex conjugate. From the analysis 
of section 4.4 the zero modes can have can have a 1/{z — Zq) pole at z = Zq and a 1/{z — Zq) 
pole ai z = Zq . 



The most general zero mode, with such singularities, satisfying ( |6.1|) is given by a 
simple generalization of 



y = -b ^'[^ -b* ^'i^ - -ib + bn^, for (6.4) 

where the h and the h are complex numbers. The numerators of (^]5) have zeros at 
z = Zq and z = Zq (see equations (|2.6|) ). The denominators of the holomorphic and 
anti-holomorphic parts have zeros at respectively z = z^ and z = Zq = z^ (see equations 
( p.l3|)(|2.14|) ). So in total the zero mode ( |6^ ) has the correct pole structure. Choosing the 
standard positions for three of the incoming and outgoing strings 2:1 = 0, 2:2 = 1, 2:4 = 00 
and = A we can write the zero mode as the simpler expression 

r = -(6 + 5-)f + (f -f)f6^ + S-^). (6.5) 

iV4 ViVi iy4^/\z — ZQ z — Zq j 

We could, of course, have written down the overall form of this mode directly from the 
constraints on its singularities specified above. The advantage, of the more complicated 
form ( |6.4| ) is that we can read off immediately its behaviour asymptotically far down the 
strings which we indicate on the r.h.s of ( |6.4| ). 

The holomorphic term in ( |6.5| ) corresponds to translation of the branch point at Zq 
and the anti-holomorphic term to the translation of the branch point Zq . To see this more 
explicitly note that in the neighbourhood of one of the zq^s we can write 

w = wq -\ — a{z — zof' -\- ■ ■ ■ and X = X{zq) + (3{z — zq) + ■ ■ ■ (6.6) 
2 



Inverting the first equation to write z — zq = ^Jljotiw — wq) and substituting into the 
second we arrive at an expression for X in terms of w. We can then differentiate with 
respect to the branch point position to find 

a 1 



9^0^ = - + (6-7) 

OL Z — Zf^ 
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(6.9) 



with a and /? given by 

a = dlw\z=z^ and ^ = d^X\z^^^. (6.8) 
Evaluating a and /3 explicitly one finds 

The final outcome of this analysis being that one can identify the shift in the branch points 
bw^ and with the parameters h and h via 

b = Swq and b = Swq . (6.10) 

Out of the four parameter space two parameters translate the whole classical solution 
in the r and a directions, and two change the relative r and a separations of the branch 
points. 

Finally, in addition to the non-trivial bosonic zero modes constructed above, there are 
eight constant modes corresponding to globally translating the world sheet in the eight 
dimensional transverse target space. 

6.3. Fermion zero modes 

From section 4.5 the pole structures allowed for the fermion zero modes are a 1/ {z—Zq) 
pole for 9j a l/{z — Zq) pole for 6*, a l/(^ — Zq) pole for 9 and a l/{z — Zq) for 6. 

The construction of the fermion zero modes is thus virtually identical to that for the 
bosons; they are given by 

0p^(P±-(PL-P±)^^)0P ^ PL0P (6.12) 
0P^(P±-(PL-P±)^l^)0P ^ PL0P^ (6.13) 

^N^ nJz-z^' Ni ' ^ ' 

~ _ ~p 

where 6^,6^,6^ and 6 are constant two component complex spinor variables. The simpler 
expressions are written for the standard choice zi — 0,Z2 — 1,2:4 = 00 and z^ — A. 
Asymptotically far down a string {z —>■ Zi) the fermion zero modes tend to the constant 
values indicated on the r.h.s. of the above expressions. As we will see in section 8 these 
asymptotic values determine the kinematic structure of the scattering amplitudes. 

In addition to the non-trivial fermion zero modes constructed above, there are eight, 
globally constant fermion modes 6^ ,9^ ,9^ and 9 which take the same constant value on 
each of the four strings. 
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7. Fermion zero modes and graviton wavefunctions 

In this section we will focus on the fermionic zero modes for a single string. This is 
equivalent to studying the zero modes of the previous section asymptotically far down one 
of the strings. We will construct ground state wave functions from these modes. 

The decomposition of the fermionic variables ( |3.3[ ) into coordinates and momenta 
explicitly breaks the spin{8) invariance down to U{1) U{1) SU{2) SU (2). Using the 
transformation properties of the 9 under SO [8) rotations however one can build up from 
combinations of the 9 the components of the 8s, 8c and 8„ representations of SO{8). We 
will restrict our attention here to the vector representation 8„ and will construct one set 
of vector components from 9^ and 9^ for the left moving sector and another set of vector 
components from 9^ and 9^ for the right moving sector. Tensored together they will then 
form the ground state quantum numbers for the 5*0(8) polarizations of the graviton G^"^, 
antisymmetric tensor field B^"^ and the dilaton $. Specifically, with polarization tensor 
e^"^ the wavefunction reads 

^'(e^^, 9, 9, 9, 9) = ^\9, 9) ^^{9, 9) , (7.1) 

with \1/ and \E' denoting respectively the left and right moving parts of the wavefunction We 
will just focus on the left moving sector, since (before integration over the moduli) the left 
and right sectors are completely decoupled and the analysis for the right movers is identical 
to that for the left movers. As a book keeping device we will contract the left moving 
components with a polarization vector ^. Using the 5*0(8) transformation properties of 9^ 
and 9^ enumerated in section 3 and the appendix the 5*0(8) vector separates into three 
parts corresponding to the subspaces (X^,X^), {X^^X"^) and {X^ , X^, X'^ , X^) 

~ ~ - - (7.2) 

= ^'(12) (e, r , 0, 9) + Vl>(34) (e, r , 0, 9) + VE'(5678) (C , ^, ^), 

where the indices indicate the corresponding subspaces. The \E'(i2), ^(34) and \I/(5678) are 
given by 

^(i2)(^,r,^,^)= ^p+{eo'+^9') ^ e=-ie±^e) 

^'(34)(l,r,^,^)= f+(p+)'r^'^' with ' 2 (73) 

^i,e7s){r.0,9)= p+rOa'^9, ^'^ 2^^ ^' 

28 



The standard definitions for 9'^ and 6a^6 are given in tlie appendix. Note tliat to obtain 
tfie correct U{1) cfiarges for the components of the polarization vector ^ the wavefunction 
\E' must be assigned the U{1) charges (0, 1) with respect to rotations in respectively the 
X^,X'^ and X^jX'^ planes. Each 9 is weighted with a factor of \/p^, where is mea- 
sured incoming to the process. In other words for the two strings 1 and 2 is replaced 
respectively by A'^i and N2 whereas for the two outgoing strings 3 and 4 p'^ is replaced by 
— A^3 and — A^4. The powers of p"*" are needed to reproduce Lorentz invariant results.^ 

Finally note that in light cone gauge an arbitrary physical polarization vector ^ can 
be constructed by specifying only its components in the transverse space. The equations 
of motion impose the constraint 

k.^ = 0. (7.4) 

There is a remaining gauge degree of freedom which leaves the physical polarization con- 
dition ( [7.4|) unchanged and corresponds to adding a multiple of momenta k to By gauge 
choice one can set the -|- component of the polarization to zero. We thus have the physical 
polarization given by 

^ = {^+,r,td,t,r), (7.5) 

with 

^+ = and ^- = }-{^p*+Cp). (7.6) 

Likewise an arbitrary physical polarization tensor e can be specified by its transverse 
components. 



8. The four graviton scattering amplitude 

The calculation involves integrating over all modes. The non zero modes lead to 
determinants, whereas the zero modes must be soaked up into the wavefunctions, or, in 
the case of the gauge field zero mode, will be integrated over a finite interval. In this 
section we study in detail the integration over the zero modes. At the end we return to 
the fluctuation determinant and combine the results to obtain the four graviton scattering 
amplitude. 

We deduce the necessary factors, by observing that the resulting scattering amplitude is 
Lorentz invariant. A more rigorous way of determining them, however, would be to construct the 
5*0(8) rotation generators out of the 6 and A and verify that, with the relative factors of ( [7. 3D , 
^'^ is indeed an £10(8) vector. We do not pursue this here. 
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8.1. Integration over hosonic zero modes 

In section 6.2 we constructed four non-trivial zero modes for the bosonic fields. In 
addition to these modes there are eight constant bosonic zero modes corresponding to 
shifting the whole string configuration through the eight dimensional transverse space. 
Integration over these zero modes leads to transverse momentum conservation. 

The non-trivial bosonic zero modes of (|6.4|)(|6.5|) correspond to translations of the 
branchpoints. Two of the possible four modes, b = b = 6t and b = b = ia correspond to 
translation of the whole classical solution in the r and a directions respectively. Integration 
over the first of these modes leads to conservation of the light cone energy, p~ . Integration 
over the second leads to invariance under shift of a 

The two others correspond to the relative displacement of the branch points. The 
integral over these modes is not suppressed by a gaussian term in the action. We thus 
have to integrate over them carefully. 

Let us denote by Wr = Tr + iar the point in moduli space of the four string scattering 
diagram. It specifies the relative displacement of the branchpoints with respect to each 
other. The point = is the saddle point world sheet. For small Wr we can identify 
Wr with the parameters b and b of ( |6.10|) i.e. Wr = b = —b. Denote Xc{w,Wr) the 
classical matrix string field configuration corresponding to the point Wr. Varying Wr from 
Wr to Wr + 6wr takes us from one classical configuration to another. The change of fields 
corresponding to this deformation are, by the equations of motion, zero modes of the 
quadratic operator; 

Yr{w,Wr) = dr^Xc{w,Wr) and Yc,{w,Wr) = d^^Xc{w,Wr). (8.1) 

However, zero modes which tend to a non-zero constant value asymptotically far down a 
string will not decouple from the action. The background field configuration ( |2.9| ) is linear 
in r and in integrating the kinetic term {d{Xc + Y))'^ by parts one picks up a boundary 
term, which couples the zero mode linearly to the classical background. For the saddle 
point world sheet configuration the zero mode corresponding to relative displacement of 
the branch points tends to zero asymptotically down the strings, and hence does decouple. 
This is another way of seeing that the background is a stationary point of the action. It 
also indicates that one needs to switch from integrating over this mode to integrating over 
the modular parameter Wr. 

A general field configuration can be decomposed in two ways. We can set Wr = 0, 
i.e. start with the saddle point classical configuration, and expand using a complete set 
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of modes about this background. Alternatively we can keep Wr as a free complex variable 
and expand in the set of modes excluding the zero modes IV and of (|8.1D . We thus 
have the identity 

X{W) =X,{W, 0)+J2 VnYniw, 0) 

V (8-2) 

n^Ta 

We know the integration measure for the first decomposition. It is determined by the 
condition. 

j (rfy)e-^<^l^> = 1 with < Y\Y' >=^ j (fwTi[Y{w)Y'{w)+Y'{w)Y{w)]. 

(8.3) 

Substituting in the first decomposition of ( ^.2|) we find the zero mode measure 



j idY)o = j l[dy,\det < Y^\Yr >\K 



(8.4) 



To obtain the integration measure for the second decomposition we need to calculate 
the jacobian for passing from integration variables x^-^Xfj and Xn to integration variables 
Triar and x^, (n 7^ r, a). 

From ( |8.2| ) we have the identity 

Y,yMw)=X{w,Wr)-X{w,Q)+ VnYniw), (8.5) 

where the sums are over all modes. Taking inner products we can relate the two sets of 
variables 

_ < Yr\X{w,Wr) -X{W,0) > 



< Yr\Yr > 
_ < Y^\X{w,Wr) - X{w,0) > 

< Y^\Y^ > 

<Yn\X{w,Wr)-X{w,0)> , _ „ . 

Vn = + Vn for n^T.a 



(8.6) 



The jacobian for passing from one set of variables to the other is seen to be given by the 
finite dimensional determinant 



J = det ( l^^y^ f^y^ ) ~ 1, (8.7) 
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where we have used (|8.5| ) (pJ] ) to see that, to lowest order, J = 1. Corrections to this 
approximation form a power series in r and a. In the domain of vahdity of the one loop 
approximation, where the momenta p are very large, these terms can be dropped since the 
integrals over r and a are entirely dominated by their quadratic corrections to the minimal 
action. The conclusion is that the integration measure for the collective coordinates r and 
a and for the other zero modes is exactly as for ( |8.4| ) except with dyr replaced by dTr and 
dija- replaced by dar ■ 



J {dY)Q = j dTrdar W dyq\dei < Yq\Yr >|^ 



By the choice of the decomposition (|8.2|) the only dependence of the amplitude on the 
relative displacements and is through the classical action. Expanding the classical 
action up to quadratic order in Tr and ar the integrations over and cr^ can then be 
performed. The integrals have to be evaluated by the method of steepest descent (see [|T2| ) 
leading to a factor of z, the final result being 

dTrdar e^^""") = i — e-i(«i°g^+"°g*+^i°g'^), (8.9) 
sut ' ^ ' 



with s,t and u the Lorentz invariants defined in (|2.30|)(|2.31D . The non-Lorentz invariant 
factor c is given, up to numerical factors, by 



c 



s^Nl\z+-z^\\ (8.10) 



where Zq and Zq are the two branch point positions in the complex z plane. They are 
determined by the quadratic equation ( p.6| ) written with the standard choice = 0, ^2 = 1, 
^3 = A, ^4 = oo, with A set to the value corresponding to a stationary point of the classical 
action (see for details). 

The non-Lorentz invariant factor c might at first appear problematic. However pre- 
cisely this factor occurs in a light-cone gauge calculation for the same amplitude in super- 
string theory or bosonic string theory where it is well known that light-cone gauge produces 
Lorentz invariant results. It is thus instructive to see how Lorentz invariance is restored 
in the case of light cone string theory. Described in detail in chapter 11 (see in particular 
appendix 11 A) of Green Schwarz Witten |T^ is a careful evaluation of the fiuctuation de- 



terminant (involving a regularization of the branch points) using the trace anomaly. The 
result being that the determinant contains precisely the right Jacobean factor to convert 
the integral J dTrdar into the integral j d^X where A is the standard complex moduli for 
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the four string scattering process already introduced in the previous paragraph. In terms 
of the complex moduli A the classical action (|2.11 ) reads 



S = ki.ks log |A| + k2.h log |1 - A|. (8.11) 

Integration over A clearly leads to a Lorentz invariant result. 

In the light of the above discussion it seems plausible that the matrix string determi- 
nant for the amplitude studied in this paper could likewise generate the same Jacobean 
factor. To show this explicitly, however, would involve finding some generalization of the 
the trace anomaly techniques applicable to the non-abelian determinant of (4J.). We do 
not address this question more fully here. 

8.2. Integration over gauge field and ghost zero modes 



There are three nontrivial gauge boson modes (|6.2| ) , corresponding to the fact that 
a generic gauge field zero mode takes on four different values, ff-, down the four strings, 
subject to the constraint J2i ^i^^i = 0- 

As for the bosonic zero modes of the previous section their integration measure includes 
a determinant of their inner products. 

{dV)o = \\ dai\del<V^\Vj>\^ 5{ai + a2-ar^-a^), (8.12) 

J io 



where the determinant is a three by three determinant and can be calculated by choosing 
any linearly independent basis for the Vi. One such basis would be to say that Vi {i = 1, 2, 3) 
is the mode which tends to 1 down the zth string, tends to ±1 down the fourth string and 
tends to zero down the others. The integrals only go up to Oi = gs since there is a globally 
well defined gauge transformation that identifies = gs with = 0. Specifically for the 
string i the globally defined unitary matrix is U = e~'^"/^\ 

Integrating over the gauge field zero modes gives rise to the contribution 

{2ngsf |det < Vi\Vj >\K (8.13) 

There is a single ghost field zero mode, the constant mode. It comes from the fact 
that there is a globally constant U{1) gauge transformation. To understand how to deal 
with this mode, let us look carefully at the Fadeev Popov representation of the gauge fixing 
determinant. We insert into the functional integral the factor of 1 written as 

1 = A l{dU)5'{f-D,V^^^), (8.14) 
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where / (dU) is the Haar measure for the integral over the unitary group and Va^^ is the 
gauge transformation of Vq, with respect to U. The prime on the functional delta function 
means that the constant mode is excluded. By the standard argument, once this factor is 
inserted into the functional integral, a gauge transformation allows one to factor out the 
volume factor J [dU) which one can then divide out, and one is left with just the Jacobian 
factor A. / in (|8.14[ ) is some arbitrary function which can be integrated over with gaussian 
weight to induce the standard gauge fixing term 

To evaluate A we expand the unitary integral about the point where the argument of 
the delta function is zero. We then have, under a gauge transformation. 
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,U^D„U = D„A with U 



e 33 



(8.15) 



We now write A, the unitary integral / {dU) and the functional delta function in terms of 
modes. 



H'^) = ^^K^nio-) with A„ = 

n 

J{dU)= j{dA)= jl[dXr,\Ar 

5'{A~B) = J]5(a„-M-^. 



< An|A > 

< A„IA„, > 



(8.16) 



We have used the flat space definition of the measure (|8.4|) for the functional integral over 
A. The normalization factors |A„| are given by 



|A„| =< AnlAn >2 . 



i.l7) 



In the expression for the delta function of ( p.l6| ) A and B are arbitrary functions which 
we have expanded in terms of the modes A^, i.e. A{a) = J2n'^nAn{o') and similarly for 
B. Note that there is nothing special about the decomposition into the modes A„, any 
(orthogonal) set of modes can be used. 

Substituting in the expressions ( |8.16| ) into ( |8.14| ) we find 



1 = A 



rfAoiAol n 



< A^\An > 

< A„,|D2|A >■ 



iA8) 



The integral over Aq is decoupled from the functional delta function and hence is not 
localized. It corresponds to a global U{1) rotation and from (|8.15| ) we see that it must 
be integrated over the finite interval [0, 2ngs]. The product (determinant) in (|8.18| ) can be 
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represented in the standard way by introducing a ghost action and integrating over aU 
ghost modes except the zero mode. We thus have 

1 



2T:g,\ < Ao|Ao > |^ 



{dcdcYe 



^.19) 



where the prime indicates that the integral excludes the zero mode. The combined contri- 
bution from the ghost and gauge zero modes is thus 



det < ViWi > 



(8.20) 



< Ao|Ao > 

The factor of g"^ is the correct Qs dependence expected from string theory with each branch 
point contributing a factor of Qs- 

As was originally argued in the zero mode structure of the abelian gauge and ghost 
fields will always generate the correct Qs weight for the topology of the string diagram. 

8. 3. Integration over fermion zero modes 

In section 6 we constructed eight non-trivial fermion zero modes which tended 
to constant but different values asymptotically far down each of the four strings 
( (|6. 1 1|) ( |6. 12| ) ( |6. 13| ) ( |6. 14] ) ) . In addition there are eight constant fermion zero modes taking 
the same value on each string, 9^ . We thus see that the fermion zero mode down the ith 
string can be written as 



Or 



1 



p,6P + NiO^ 



Z.N 



(8.21) 



\p*9 +N,9 ) 



We now integrate over the fermion zero modes including a graviton wavefunction for each 
of the four strings. Since the left and right moving sectors of the fermionic integrals are 
identical we only need focus on one of them. We have the following zero mode integral for 
the left moving sector 



M / d^d n^'(^-^-^^) ^^^^ 



dH=j dH^ d^QP dH"" d^d"" 
A/'=|det < rir >r^ 



.22) 



M is the normalization for the zero mode integral. It follows from the fermionic equivalent 
of ( |8.3|) . The wavefunctions are given in ( |7.2| ) and the zero modes Qi and Qi are given in 



terms of 6'p fi'P fi^ , and 6'^ through equation (|8.21| ). The fermionic integrals will pick out 
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from the product of wavefunctions precisely two 6'^ , two 6^ ^ two 6^ and two 6^ . Using the 
fact that the wavefunctions break into three subspaces (|7.3D the fermionic integrals pick 
out of the product of wavefunctions six different types of contribution. In other words, 
under the fermionic integrals we have the identity 

4 4 4 4 

n ^^(6, = n ^a2) + n ^(34) + n ^(5678) 

i=l i=l 1=1 i=l 

+ E (^(12)^(12)^(34)^(34)+ ^(12) ^('l2)*(5678)^(5678) 

+ *(34)*(34)^(5678)^(5678)) 

(8.23) 

Each of these six contributions has to be evaluated separately. The calculations are lengthy 
but straightforward. The appendix lists the identities necessary to derive the results and 
presents the first calculation in detail. It also gives essential intermediate results to aid 
verification of the five other calculations. Although the calculations are not particularly 
enlightening in themselves, the reader is nevertheless urged to turn to the appendix to 
convince him(her)self that it is highly non-trivial that the fermion zero mode integrations 
produce the ten dimensionally Lorentz invariant results enumerated below. In this section 
we will simply state the results, referring the reader to the appendix for all details of the 
calculations. 

We start by considering the case in which all polarization vectors are in the 
jX"^ space. Performing the fermionic integrations we find, after several pages of 
algebra, (see appendix) that we can express the final result in terms of Lorentz invariant 
quatities. Specifically we find 

/ d'9l[^l,,js) = -h^(6-6)(6.e4) + st(ei.6)(6.e4) +si^(ei.e4)(6-e4)], (8.24) 
^ 1=1 

where s, t and u are the ten dimensional Lorentz invariants given in (p. 301 ) ( p73l|) . This is 
of the form expected for all four polarization vectors orthogonal to the momenta, (see for 
example ^^0) 

The next case considered is that in which all four polarization tensors lie in the X^,X'^ 
space. Integrating over the grassman variables one finds, after several pages of algebra (see 
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the appendix for essential details) that the result can also be expressed in terms of Lorentz 
invariants. The final result reads 

4 

1=1 

4 [s^iii2i*3ii + u%isi;ii + t'i^ut2i3 + t^isiiii + u^i2idii3 + s^isueM ■ 

(8.25) 

Using the identity ( |2.32| ) and the definitions for ^ and ^ of (|7.3|) it is straightforward to 
show that this is again of the form ( |8.24| ). 

Focusing next on the case with two polarizations in the space and two 

in the X^^X"^ space we find that there are six types of contribution depending on the way 
the pairs of polarizations are divided amongst the four strings. Performing the fermionic 
integrations we find after a long calculation (see the appendix for essential results) that 
the results can be expressed in the Lorentz invariant form 

I rf'^*;34)*(34)*(5678)*(5678) = -2t^i(lll2 + I2I1*) (6-^4) 

I rf'^^;34)*(34)*(5678)*(5678) = -2«t(fll3 + Isfi*) (6-^4) 

J rf'^^;34)^(34)*(5678)^(5678) = -^Sui^lQ + ^€1^2.^ 

J ci«^vE'?34)^'j34)*;5678)*(5678) = -2««(6f3 + 6^2*) (6-6 

J CJ«^*?34)*J34)*(5678)^(5678) = -Sstds^ + I4I2*) (6-6) 

J C^'^^'?34)*(34)^(5678)^(5678) = -2«t (6^ + 144* ) (6 -6) , 

where, s,t and u are the Lorentz invariants defined in (|2.31|) . This is again of the 
form ( |8.24| ). In conclusion we see that for arbitrary polarization vectors lying in the 
X^ ,X^ ,X^ ,X^ ,X''' ,X^ space, i.e. transverse to the momenta, the result can always be 
written in the form (|8.24|) . 



3j 



4 
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We still have to evaluate the fermion integrals for the three cases where two or more 
polarizations are in the X'^,X'^ plane. We begin with the case in which all four polarizations 
lie in the X^,X'^ plane. Evaluating the fermion integrals we find a relatively simple result 
(see the appendix). The result can then be shown to be equal to the Lorentz invariant 
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quantity given below. 



4 



1=1 

+ (6-fc3)(?4.fc2)(6-6) + (6-fc4)(6-fcl)(6-^4) 
+ 2t[{^2.kl){U-hmUs) + (6-^4)(6-fc2)(6-e4) 

+ {C2.k4Mi.ksM3.U) + {^3.ki){U-k2m-^2] 

+ 2u[{^i.k2){U-k3){^2.^3) + {^3.kM2-kim-U) 

+ {^i.k4M2.ksM3.U) + {^s.k4){U-kimU2] 



1.21) 



The vector products are the Lorentz invariant vector products, with ki the ten vectors 
of momenta ( |2.25|) and ^ the ten vectors of polarization which include a ^~ part given 
by equation ( |7.6D . To verify equation ( p.27| ) one substitutes into the right hand side the 
mometum ten vectors ( |2.25| ) and polarization ten vectors (|7.5|)(|7.6| ). The complicated form 
of the right hand side of ( ^.271) can then be shown to reduce down to the simple result 
given in the appendix. It is recommended to use a computer algebra program such as 
Mathematica to check this. 

The result ( |8.27| ) is precisely the four graviton kinematic factor calculated in string 
theory |^^[|T9[. It is totally symmetric under interchange of the strings and is further- 



more gauge invariant. In other words if one replaces any one polarization tensor by the 
corresponding momenta the result is zero. 

Next we treat the cases in which two polarizations lie in the X^,X'^ plane and two 
in the X^^X"^ plane. There are six different possible ways of distributing the two types of 
polarization amongst the four strings. Computing the fermionic integrals one arrives at 
some relatively simple expressions (see the appendix). It is then straightforward to show 

38 



/ ^'^^(12)^(12)^(34)^ 
/ ^'^^(12)^(12)^(34)^ 



(34) 
4 

(34) 
3 

(34) 
4 

(34) 



that they can be expressed in the Lorentz invariant form hsted below. 

j d'e<fl,,)<fl,)%,)%,) = -[tlt(ei.6) -2t(6-fc3)(6.fc4) -2z/(6.fc4)(6.fc3)J(6-e4) 

-2s(6.A;4)fe.A;2) -2t(6.A;4)(6-fc2)](6-^4) 

-[tis(a.e4) -2s(a.fc3)(e4./C2) -2tl(6.fc2)(e4./C3)](6-e3) 
-[St(6.^4) -2s(6.fc3)(^4.fcl) -2t(6.fcl)(^4.fc3)](6.6) 



3 

(34) 
2 

(34) 



These are again of the form ( fj.27| 
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FinaUy for the case with two polarizations in the X^^X"^ plane and two in the sub- 
space we find (see appendix) identical results to those of ( p.28|) with the 
wavefunctions \I/(34) replaced by \I/(5678)- 

The result of this analysis is that an arbitrary physical polarization vector, determined 
by its components in the eight dimensional transverse space, leads to the Lorentz invariant 



result of equation (|8.27|) . The restoration of the full SO {8) invariance of the theory was to 
be expected from a correct quantization of the theory, but the extension of this invariance 
to ten dimensional Lorentz invariance seems to be a very non-trivial test of matrix string 
theory. 

8.4- The four graviton scattering amplitude 

Combining the integrals over the left and right moving fermion zero modes, the bosonic 
and gauge zero mode integrals and the determinant from the non zero modes of section 4 
we arrive at the final result for the four string scattering amplitude : 

,2 

'^AtlM2M3A'4'^I'lI^2!^3l'4tl t2 ^3 ^4 ' {O.^^yj 



9 s 



j^ — 2_^J_ g-^(slogs+tlogt+-ulog-u) 

stu 



where the tensor JC 



MiA'2M3A*4 



is the coefficient of ^f'^f .^f in equation {^^l7^. The 



factor X contains the finite determinants from the zero mode normalizations along with 
the fiuctuation determinant of section 4. Up to an overall numerical factor it is given 
by 

< AolAo > det < ri^" > 



X= cj 



det < VAVi > det < YJYr > 



130) 
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where c is the non-Lorentz invariant quantity defined in ( |8.1C1|) . We have shown in section 
4 that the fluctuation determinant J does not receive any singular contributions from the 
branch points. We further argued that up to smaU corrections the bosonic and fermionic 
determinants would cancel, but we do not know how to evaluate J more precisely. We can, 
however, take some inspiration from conformal fleld theory. In CFT it is well known that 
the flnite determinants for the zero modes combine with determinants from the non-zero 
modes to give a result (in the critical dimension) that is conformally invariant. The heat 
kernel methods used to show this (see for example [^) do not generalise in any obvious 
way to the determinants calculated in this paper. It seems likely, however, that a similar 
mechanism is at work for the combination of determinants of ( |8.3C1| ).^ As already discussed 
at the end of section (8.1) it is also quite possible that the non-Lorentz invariant factor c 
would also be absorbed in a natural way into the determinant (again this is precisely what 
happens in the string theory calculation for the same amplitude). 

If this is the case then, up to an overall numerical factor, we have reproduced from 
matrix string theory the high energy limit of the string theory scattering amplitude for 
four gravitons pi [O P] . 



9. Discussion and conclusions 

In the previous section we have performed an analysis of graviton scattering. This 
formalism can however be extended in a straightforward way to wavefunctions involv- 
ing fermionic ground states for the left and/or right moving sectors. These would be 
constructed from odd numbers of ^, and would permit the construction of the complete 
multiplet of massless states of the type IIA superstring. 

9.1. Summary of the general calculational procedure 

We study scattering amplitudes which in string theory are dominated by the points in 
moduli space corresponding to saddle points of the classical action. We focused on the tree 
level contribution, but in string theory there is a whole tower of saddle point world sheets 
of higher topology that also contribute.^ The methods developed in this paper should 

^ I thank Herman Verlinde for emphasising the importance of this effect in CFT and for 
suggesting that a similar mechanism could play a role here. 

^ In the original analysis of Gross and Mende |12} | it was found that the perturbation series 
was strongly divergent with genus g world sheets giving a contribution proportional to g^^ . 
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be directly applicable to these higher genus contributions. Let us summarize the general 
procedure. 

The classical matrix string solutions corresponding to the saddle point world sheets 
are used as background field configurations. The quantum fluctuations around the back- 
grounds are treated by a one loop calculation. This consists of a zero mode part and a 
fluctuation determinant for the non-zero modes. All the essential structure of the am- 
plitude is contained in the zero mode integrations which can be calculated exactly. To 
flnd the zero modes one flrst uses the fact that around the interaction points there are 
instanton like fleld conflgurations. Locally these break the translation symmetry and part 
of the supersymmetry. Each broken symmetry has a corresponding zero mode. They are 
given by equations ( |5.13| ) ( |5rT9| ) 

V^ = F^^ and 6 = F^,V^,e. (9.1) 

where F^i, is the background field strength written in ten dimensional notation (see section 
3.2). In addition there are zero modes for the two dimensional gauge field. They can be 
written as 

= D^Xi, (9.2) 

where Xj is any of the non zero background bosonic fields. The zero mode field config- 
urations (pJ] )(|9.2| ) are finite at the interaction points. Since, however, the background 
bosonic field configurations Xj correspond to world sheets branch points, the zero modes 
( |9.1| )( |9.2D will behave, asymptotically far from the branch point, like singular string world 
sheet zero modes, with the singularity sitting at the branch point. These modes can be 
glued into global world sheet zero modes. 

The important point about the construction of these global zero modes is that they 
only depend on the existence of a finite instanton like classical solution around the branch 
point, not on the precise details of the field configuration. 

For any world sheet it should be possible to identify, using (|9.1j)(|9.2|), the allowed 
singularities of the zero modes at the branch points and hence to construct the global zero 
modes. 

Integration over the gauge field zero modes, along with the single, constant, ghost 
field zero mode, leads to the correct power of Qs for the string world sheet The global 
bosonic zero modes correspond to the collective coordinates for the positions of the branch 
points. Integration over the global fermion zero modes gives the ten dimensionally Lorentz 
invariant kinematic structure of the scattering amplitude. This we have has explicitly 
verified for the case of four graviton scattering. 
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9.2. Connection with light cone superstring calculations 

The matrix string calculations described above have a strong resemblance to the func- 
tional integral methods used in light cone superstring calculations (see the review article 
||25|| ). Below we discuss the similarities and differences between the two. 

In light cone superstring calculations the functional integral reduces down to an inte- 
gration over zero modes. There are, however, no zero modes corresponding to the trans- 
lation of the branch point and no abelian gauge field and so no gauge or ghost field zero 
modes. The inclusion of factors of gs and the integrals over the string moduli are added 
by hand. 

To perform the Euclidean functional integral the Minkowski space Majorana-Weyl 
fermions are combined into complex fermionic coordinates and their conjugate momenta. 
These can be defined by the same formula as for matrix string theory i.e. equation ( |3.3| ). 
This breaks the spin{8) symmetry down to U{1) (S> SU{4). The 9 and 9 of ( p.3| ) form a 4 
of SU{4) and the 9 and 9 form a 4. 

In the matrix string theory calculation, there is a classical background, that through 
the non-abelian nature of the theory, interacts with the fermions and further breaks the 
SU{4) symmetry. In the four string scattering process discussed in this paper the original 
spin{8) symmetry is broken down to U (1) (8> U{1) ® SU{2) (g) SU (2). Other classical back- 
grounds, for example five or six string scattering, would break the spin{8) symmetry even 
further. 

In light cone superstring calculations the fermionic coordinates 9 and momenta A have 
the conformal weights 1 and respectively. The functional integral over the fermions leads 
to a determinant and to an integration over the 9 zero modes. These are given by the 
complete set of abelian differentials defined on the string diagram Riemann surface. In 
terms of the light cone world sheet coordinate w, they consist of all the modes satisfying 

d^9 = 0. (9.3) 

We denote by 9 in ( |9.3| ) the 4 of SU (4) made from 9 and 9 of ( p73|) . There is also the 
complex conjugate equation for the 4. The modes satisfying (|9.3|) are allowed to have 
1 / y/w — wq singularities at each simple branch point wq, and have the boundary condition 
that they tend to constant values asymptotically far down the strings. The generic fermion 
zero mode thus has a singularity at each branch point. This is in contrast to the matrix 
string zero modes of section 6.3 which each had a single singularity even though there were 
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two branch points. In other words, for this process, there are twice as many fermion zero 
modes in the hght cone superstring calculation as there are in the matrix string calculation. 

This mismatch is resolved by the final element of light cone superstring calculations : 
the addition by hand, at each branch point of an operator V needed to preserve supersym- 
metry. The precise form of this operator for simple branch points (i.e. where two strings 



interact) was calculated in [^. It can be written in the form 



^ = - 10^2^/ M "^{dzX', dj, dj) ^{d,x', die, die). (9.4) 

where \I/( , , ) is as in ( [7^21) except with p+ replaced by l/(2d'^w{z)), and in the right 
moving part \E'( , , ) p+ is replaced by l/{2d'^w{z))* . w is the light cone coordinate 
expressed, through the Mandelstam mapping, in terms of the uniformization z (for tree 
level scattering see equation ( ^.5| ) ) . The bosonic fields X in the operators V are contracted 
both with the external wavefunctions and with the other V. The fermions e are the fermion 
zero modes. It is straightforward to see that the two interaction operators used for the 
four string scattering amplitude soak up the excess zero modes. 

9. 3. Conclusions and future directions 

The matrix string zero mode calculation of this paper is considerably simpler than 
that for the light cone superstring, where an interaction operator has to be introduced 
by hand at each interaction point. So far, however, the matrix string calculations have 
only been carried out about a classical background corresponding to a saddle point of the 
classical action. It is important to extend them to other backgrounds. In the context of 
the four-string scattering calculation, it is clear that there is a problem when the world 
sheet is no longer a stationary point of the action. In this case, the world sheet is not 
(anti)holomorphic around the branch points and the corresponding instanton solution will 
break all of the supersymmetry leading to extra fermion zero modes. There must be a 
simple mechanism that soaks up these extra zero modes. It is important to understand 
this as it could permit the calculations to be done at an arbitrary point in moduli space. 

The semi-classical analysis of this paper has been justified by focusing on very high 
energy scattering processes. Away from this limit we do not know how to calculate. There 
will however be an effective action describing the low energy fiuctuations, even if we do 
not at present have the technical tools to obtain it. The zero mode analysis of this paper 
only depends on the existence of a finite field configuration for the interaction points, and 
should also be applicable to this effective field theory. 
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Finally it is interesting to speculate on how the semi-classical reasoning of this paper 
could be made more rigorous. In simpler situations, semiclassical methods have recently 
been shown to be in precise agreement with exact results. Using quasi-classics, the partition 



function for matrix string theory on a torus was calculated |29] and shown to agree, on 
extrapolation to the small torus limit, with the exact results for the completely reduced 
SYM theory [^^. These latter results were obtained by mapping the original SYM theory 
to a cohomological field theory (CohFT) using the methods of |31]. Recently CohFT 



methods have been applied directly to the matrix string partition function beautifully 
confirming the semi-classical results of |5^. A key aspect of CohFT is that the functional 
integral localizes on configurations of the classical vacua. Although it is not possible to 
apply the twisting procedures used in |^ to matrix string scattering processes it is 

tempting to speculate that some kind of localization mechanism is also behind the success 
of the calculations presented in this paper. 
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Appendix 1. Gamma matrix definitions 

We choose as the basis of real spin(8) gamma matrices F^ and Majorana-Weyl fermions 

where the indices a and a which run from one to eight label the 8s and 8c representations 
of spin{8) respectively. The 7^^ (which are the transposes of the ■jI^) are given by 

7^ = l(g)l(g)l 7^ = e(g)l(g)l 

7^ = Ti (g) Ts ® e 7"^ = Ti (g) n ® e 

(1.2) 

7^ = Ti (g) e O 1 76 = T-3 (g) e O Ti 



7'' = — Ts (g 1 (g e 7^ 



Ts (g e (g Ts, 



where e = iT2 and Ti,r2 and are the pauli matrices. With this basis the gamma matrices 
r'',r^,r^ and F^ are given by 
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where all entries correspond to 8x8 blocks with the pauli matrix r2 being written as 4x4 
blocks. 

With respect to these four gamma matrices the sixteen component Majorana-Weyl 
spinors thus decompose into four blocks of four. We will be working in Euclidean space 
where ten dimensional Majorana-Weyl spinors do not exist. Furthermore, as discussed in 
section 3 each Majorana-Weyl fermionic variable is simultaneously a fermionic coordinate 
and its conjugate momenta. We thus have to combine the fermionic variables into complex 
fermionic coordinates and their distinct conjugate momenta. It is thus convenient to work 
with a complex basis of gamma matrices defined by 

J ^ {v ®v ^l)^{v^ (Siv^ ®1) with t) = -i= Q . (1.4) 

Written out explicitly they are given by 

7^ = 1®1®1 72 = -zr3®l®l 



7^ = - iT2 ®Ti®T2 7^ = iT2 ® T2 ® T2 

7^ = ^r2 Ta (g) 1 7^ = - i^i ® T3 (g) n 

7^ = — Ti (g 1 (8 r2 7^ = — ^rl (8 Ta (g) Ta , 



;i.5) 



In the basis of (|1.5| ) the four gamma matrices that appear in the background covariant 
Dirac operator are identical to those of (|1.3|) except for where T2 is replaced by — ra. 



Finally it is useful for the study of the fermion zero modes and determinant to put the Dirac 
operator into block diagonal form by mixing the two spin{8) representations. Specifically, 
writing the sixteen by sixteen F matrices as four by four blocks, we interchange the second 
and third rows and columns. This leads to the following block diagonal form for the gamma 
matrices F°,F9,F1 and F^ 



1 V (73 y \ a2 I (71 



;i.6) 



where again all entries are 8x8 blocks. This transformation also interchanges the 9^,Xa 
with ^ ,A in (2.1) of append! 
fermion action given in (|3.10| ). 



with ^ ,A in (2.1) of appendix 2 and leads to the block diagonal form for the Euclidean 
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Appendix 2. Fermionic coordinates and momenta 



Under the change of basis ( |1.4| ) the Majorana-Weyl fermions split into fermionic co- 
ordinates 9 and momenta A as follows 



^ {v®v®l)S" 



Va^/ 



and ^{v®v® l)S^ 



A 

h 



(2.1) 



where the indices A, A take the values 1, 2 with the only non-zero anticommutators being 



feA^}=5| 



and 



f2.2) 



where we have not explicitly included the spatial delta functions. The transformation of 
the fermion coordinates and momenta under rotations in the X^,X'^ and X^^X"^ planes 
and in the four dimensional space and X^ can be read off from the explicit form 

for the rotation generators in the basis ( |1.5| ). Specifically under an infinitesimal rotation 
described by w^"^ we have 



58'^ =-w''^iiS' 



1 



5^" 



with 



' n h 



7l/=^((7')V-(7')V) 



(2.3) 



We can thus draw up the following table for the action of 7^"^ on 6'"^,6'^,A^ and A 



TT 

34 
56 
57 
58 
67 
68 
78 







0-A 


A^ 


A^ 


Tg ® 1 


® i 


i 


—i 


—i 


i 


1 ® T3 


(g) i 


i 


i 


—i 


—i 


Ts 1 ® 


iri 


iri 


—iri 


—iri 




Ts (g) Ta (g) 


iT2 


ir2 


-iT2 


ir2 


-iT2 


Ts (g) 1 g 


) ITS 


irs 


-iTs 


-irs 


irs 


1 ® Ts g 






ITS 


-irs 




10 1®- 


-iT2 


-iT2 


-iT2 


-iT2 


-iT2 


1 ® r3 g 




iTi 


ZTi 


-ITi 





(2.4) 



There is an almost identical table for the action of 7^'^ on 6^, O^^^a and A . The only 
difference from the above table being that all entries of the first line change sign. We can 
read of directly from (|2.4| ) the U{1) charges for rotations in the and X'^,X^ planes. 
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Further more we can read off the effect of infinitesimal SO{4) rotations in the X^^X^^X'^ 



and X subspace. Using the notations of [p3[ for four four dimensional spinors we have 

(2.5) 



4 ^ 'A 4 
where (cr^^)^^ and (o""^")^-^ are defined by 

f ^mn\A ^I'^'m^n- ^n-m^ ( B "^/^"i^n z^ri f;^ 

(c^ ) B = 2^ (7 -a a ) and (a )^ - o" a j. (2.6) 

with a"^ and a""' defined through the Pauli matrices 

(a-)^^ = (z,ri,r2,r3) and (a"^)^^ = (-z, n, r2, rg). (2.7) 
^. i. Identities for two component 4d spinors 



We use the conventions of modified by a factor of i for the zeroth component 
since the space X^,X^,X'^,X^ has Euclidean metric. 

There are two pairs of two component spinors, 9^ and 9g with A,B — 1,2. Indices 
are raised and lowered using the antisymmetric tensor e : 

9^ = e^''9B, eA = eBA0'', ^^ = 6^%, = e^^^^, (2.8) 
where the e tensor is definied to be 

AB , AS -.1 / 1 



CAB = e"'^ = -e^B = -e with eAB=(^_-^ qJ- (2.9) 
We define 9"^ and by 

9^ = 9^9 A and = (2.10) 

Integrals over 9 and ^ are defined by 

d'^9 = J d9id92 and J = J d92d9i- (2.11) 
We define the four matrices {a'^)'^^ with m = 5, 6, 7, 8 and their barred partners by 

a^ = (z,r), and a"^ = (-z,r). (2.12) 
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Finally we contract the cr"^ with 9 and 9 to produce a four vector in the space 
and : 

9(t'^9 = 9a{<t'^)^^9b. (2.13) 

Using the above definitions we can derive identities for calculating the fermionic inte- 
grals of section 8. Below we list the complete set of identities needed. 

e^^ecB=5c and e^^e^S = (2-14) 
Q2^_^ABQ^g^ and 9^ = -e^^9A9B, (2.15) 
9^ \ { 2 \ , f P \ f 2 



and 



J 919.) = [-is) J = J ' ^'-''^ 

i^nsA = {'j'^f^ecAe^^, (2.17) 

(^m)^B(-n)^^ + (a")^^(a-)^c = 25"-5^, (2.18) 

Tr[(7™^"] = 25""^. (2.19) 



Appendix 3. Fermion zero mode integrals 

Using the identities of the previous section the calculations of the fermion integrals 
of section 8, are straightforward but somewhat tedious. In this section we present one 
calculation in detail. The other calculations can be evaluated in a similar manner and so 
for them we just give some essential intermediate results. 

For all four polarizations lying in the X^^X^^X"^ space we have 

/4 j-mi pm2 c'n^z C'^i 

n (5678) NiN2NsN4 ' ^ ^ 

where 

/• _ _ 4 

jrmm^msm^^ / d'^9Pd^9Pd'^9^ d'^9^ Y[{p*9P + Ni9^){a'^'){pi9P + Ni9^). (3.2) 

i=l 

The integrals over the 9 only gives a non-zero result for the terms from the product in which 
there are two 9^, two 9^, two 9^ and two 9^ . In other words each non-zero contribution 
will consist of two factors of transverse momenta p, two factors of transverse momenta 

48 



p* and four factors of p'^ momenta A^. We thus see that there are three different types 
of integrals to be calculated according to how the four p's and four A^'s are distributed 
amongst the four wavefunctions. Specifically we have 

J = Ti+T2+X3 (3.3) 

where Xi consists of 6 (=4!/(2!2!)) contributions of the form \pi\^\pj\'^ N'^Nf , I2 consists 
of 24 (=4!) contributions of the form \pi\'^p*pkNjNkNf and T3 consists of 6 (=4!/(2!2!)) 
contributions of the form p*p*pkPiNiNjNkNi. For these three different types of term the 
9 integrals lead to the following contributions 



\p,\^\Pj\^NiN^ j (fe 9Pa'^'9P9Pa'^^9P9^a'^''9^9^a'^'9^ 
\p,\^p*PkNjNkNf j (f9 9Pa'^^9P9Pa'^^9^9^a'^''9P9^a'^'9^ 

= -2\p,\'^p*PkNjNkN^{6'^''^^d'^^'^' ^^m^m^^m.mi _ ^m^rui ^m^m^^ 

P*P*PkPiN,NjNkNi J (f9 9Pa'^^9^9Pa'^^9^9^a'^''9P9^a'^'9P 

= 4.p*p]pkPiN,NjNkNi5^^^^5'^^^K 

(3.4) 

Adding together all the different contributions of a particular type and using the identities 
(^),(|]^) and (|OTD we find 



Ti = 4|^iViiV2iV3iV4 [27ViiV2iV3iV4^ + (A^l^l + NInI)B + [N^NI + NInI)C] , (3.5) 



.2 ^2 



X2 = 4^iViiV2iV3iV4 [j- ((ivf + Nl)N^N^ + (iV| + iV|)iViiV2) {-A + B + C) 

-q^{NiN^ + N2N^){A-B + C) (3-6) 
- ^^(iViiVg + N2N^){A + B-C)], 

and 

4 

T3 = 4iViiV2iV3iV4 [q'A + 2^NiN2N^N^{-A + B + C)\. (3.7) 
where we have defined the tensors A^B and C by 
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Adding these three results together to obtain X ( |3.3|) , substituting into (|3.1| ) and using the 
definitions ( p.31| ) we obtain the result ( |8.24|) . 



For all four polarizations lying in the X'^ ^X^ space we have 

4 



(34) ~ 2 /v-2^ -^(-*-' ^) ^ other permutations, (3-9) 

iVq iVyl 



i=l - 



where 



I{iJ;kJ) = J n (p:;^^ + iV^^^)2(p^^P + iV^^^)2. (3.10) 

m=k,l 

Evaluating the fermionic integrals (there are five diflPerent types of contribution) we find 

J(^, j; /c, = leOpfcl'iVf - 2\pk\\pkP\ +plpi)NkNf+pl{plfNlNf] + (k ^ I) 

(3.11) 

+ 64\pkf\pi\^NiN^. 

Using the identities ( |2. 2 8| ), ( 12^291 ) and ( |2.31j ) to evaluate the above result for the six possible 
permutations of (|3.9|) we obtain the result ( |8.25|) . 

For two polarizations lying in the X^^X"^ plane and two lying in the 
space we have 

/ ^^^^|34)^;34)^(5678)^(5678) = ^f^^^'^^ li^jlkj) + C.C. (3.12) 

where 

X(z;j|A;,05""'="'' = 

e J d^e (p^ep + N^e^fipjep + Nje^)^ H {p:;^ep + Nr,e^){a"'-){pjp + nJ^), 

n=k,l 

(3.13) 

where the sign factor e is +1 if /c, Z = 1, 2 or 3, 4 and —1 otherwise. This factor comes from 
the sign of the p"*" momenta in the definitions of the wavefunctions (|7.3|) . Integrating over 
the 9 (there are ten different types of contribution) leads to the result 



X{i;j\k,l)=-8e 



\Pj\'NiN^ + Nf\pk\^\pi\' + \p,f\pk\^NfN^ + \p,\'\pi\^NfNi 

Ni 

(3.14) 



+ \Pj?{PkP!+plPi)NfNkNi + {pjplpt + {p]fpm)N]NkNi 



- \Pj\^{PjPl+P*Pk)NjNkNf - \pj\\p,p;+p*pi)N,NiNi 
- \Pk\HPjPi+P*jPi)NfNi - \pi\^iPjPl+p*Pk)NfNk 
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Using the identities ( |2.28| ) , ( |2.29| ) and ( |2.31| ) to evaluate the above result for the six possible 
combinations of pairs of polarizations leads to the results listed in ( |8.26|) . 



For all four polarizations lying in the X'^,X'^ plane we have to evaluate the integral 

J <fe n vp|i2) = J(l, 2; 3, 4) +5 other permutations, (3.15) 

1=1 



where 



i{i,3-.Ki) = j (fe J] X[{piep + Nr,e''f{pmep + Nr^e''f. (3.16) 

m=i,j n=k,l 

Computing the fermionic integrals and using the identities ( |2.28| ), ( |2.29| ) and ( |2.31J ) leads 
to the following expressions for the 2{i,j;k,l) : 

J(l,2;3,4) =16N\plfpi 

J(l,3;2,4) =16[6^(NiN2NsN,f + (pDVsNI NI + pUp^f N! 



'^^NiN2NsN4{pIp3N2Ns+PipINiN^)] (3.17) 



2 



J(l,4;2,3) =16[6^{N^N2NsN^)^ + (pDYsNI Ni + pUp*,)^ Ni 



4 



Ar4 



2 

+ ^j^NiN2NsN4{pIp3N2N^ + piplNiNs)] . 

The other three I{i,j;k,l) can be obtained by taking complex conjugates of the above 
results. 

Substituting in the ten vectors for the momenta and polarizations into the Lorentz 
invariant expression (|8.27D one can show that the expression ( |8.27| ) is equal ( |3.15| ) ( |3T7| ). 
It is crucial, to obtain this equality, to include in the polarization vectors the non-zero ^~ 
part (see equation (|7.6|)). 



Turning to the case with two polarizations in the plane and two in the X^^X"^ 

plane we have 

rf8^vl/Ji2)Vl>|,2)VE'f34)Vl>j34) = .^jJL±I{r,j\k;l) + 3 other permutations, (3.18) 



where 



Iii;j\k-J) = -e I (fe n \[{p*meP + Nme'')\pJ^^Nj'')\ (3.19) 

^ ™ -7 ■ 1 



m—j,l n—i,l 
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where again e is a sign factor which is equal to +1 if z,j = 1,2 or 3,4 and —1 otherwise. 
It comes from the sign of the momenta in the definitions of the wavefunctions ( [7.3|) . 
Evaluating the fermionic integrals and using the identities (|2.28|) , ( p.29|) and (|2.31|) one 
finds that 

I{i;j\k;l)=I{i;j\l;k) and I{i;j\k;l)=r{j;i\k;l), (3.20) 

with the 6 (= 4!/(2!)^) different ways of distributing the two types of polarization amongst 
the four wavefunctions being given by 

4 2 

J(l; 2|3; 4) =16 [-^iViiV2(iV| - ANsN^ + iV|) + 2^{Ns - N^){Nipipl - N2PIP3) 



Ar2 
N2 



J(l;3 
J(l;4 
J(2;3 
J(2;4 
J(3;4 



2; 4) =16[^N,Ns + 2p^plps + 



2; 3) =16[j^NM - 2pViP3 + -^T^iPlfpl] 

.4 



iV2 



1;4) =16[^N2Ns - 2pYiPs + J^^iP^Vs] 



1; 3) =16[^iV2iV4 + 2pYiP3 + J^^iPD^Pl] 

4 2 
1; 2) =16[-^NsN4{Nf - AN^N^ + iV|) - 2^{N^ - N^m^pipl - Nsplps) 



(3.21) 

By substituting in the ten vectors for the momenta and polarizations into the Lorentz 
invariant expressions ( p.28|) one can show that the expressions of ( p.28|) are equal to those 
of ( |3.18| ) ( 13721] ). Again it is crucial to include in the polarization vectors the non-zero ^~ 
part (see equation (|7.6| )) to obtain agreement. 



Finally we turn to the case with two polarizations in the X^^X"^ plane and two in the 
X^,X^,X\X^ subspace. For this case the fermionic integrals read 

/'''»*;i2)*?i2,*f5e™,*',.6™, = ||^I(i;j) + C.C., (3.22) 

where 

j(i;j)5""='"' =- / (fe {p*ep + NjO^fip^ep + N.e^f 

_ _ (3 23) 

l[{p*JP + Nr,9^){a^-){pJP + Nj^). 

n=k,l 
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Evaluating the integrals the results are found to be identical to those of ( p.21|) , i.e. we 
have 

l{i;j) = h:{i-j\k-l), (3.24) 
where the X(i; j|/c; /) (which have /c, / 7^ z, j) are given by the expressions of ( p.21|) . 
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